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COMPARISONS OF TURBULENT -BOUNDARY-LAYER MEASUREMENTS 
AT MACH NUMBER 19.5 WITH THEORY AND AN 
ASSESSMENT OF PROBE ERRORS 

By Ivan E. Beckwith, William D. Harvey, and Frank L. Clark 
Langley Research Center 

SUMMARY 

Measurements of total temperature and pitot pressure have been made across the 
boundary layer on the wall of an axisymmetric contoured nozzle. Data were obtained at 
one station near the exit of the nozzle where the nominal stream Mach number was 19.5. 
The test gas was nitrogen and the ratio of wall temperature to settling -chamber tempera- 
ture was about 0.17. The Reynolds numbers based on momentum thickness varied from 
about 3500 to 4800. The corresponding boundary -layer thicknesses were 11.2 to 8.6 centi- 
meters (4.4 and 3.4 in.). 

The relatively low Reynolds numbers and high temperatures in the present tests 
necessitated large corrections to pressure and temperature measurements. A bare -wire 
thermocouple probe was used to obtain the total -temperature data which have been cor- 
rected for heat conduction and radiation losses. The pitot -pressure measurements have 
been corrected for viscous and rarefaction effects based on data obtained during the pres- 
ent and previous investigations. Static -pressure measurements in the boundary layer 
were corrected for viscous interaction effects. 

Comparisons have been made of the resulting Mach number, velocity, and tempera- 
ture profiles through the boundary layer with previous data and with predictions from 
finite -difference procedures for calculating compressible turbulent boundary layers. The 
agreement between predictions and the present data indicated that turbulent Prandtl num- 
bers may be considerably less than unity and that the boundary layer in the throat region 
of the nozzle may have been laminar for the present test conditions. 

Values for the exponent in the power-law velocity profile and estimated values of 
viscous sublayer thicknesses from the present data are compared with values from pre- 
vious data. Estimates of skin friction and heat transfer from the slopes of the profiles 
at the wall have been made and compared with predictions. Static-pressure measurements 
across the boundary layer indicated that the pressure increases as the wall is approached. 
Estimates of root-mean-square velocity fluctuations through the boundary layer from 
a simple mixing -length relation are compared with values obtained from the normal 
momentum equation and the static -pressure variation across the boundary layer. 


INTRODUCTION 


There is increasing interest in cold-wall turbulent -boundary -layer flows at hyper- 
sonic speeds because of their occurrence on high performance vehicles. Theories avail- 
able for calculation of turbulent boundary layers depend on empirical relations between 
the mean and fluctuating flow. For hypersonic flow these relations may be altered by the 
effects of density and pressure fluctuations. Advancements in the theoretical treatment 
of incompressible turbulent boundary layers were reported in reference 1. Several 
numerical methods, developed primarily for calculation of incompressible turbulent 
boundary layers, have been extended to compressible flows (see, for example, refs. 2 
to 4). A recent method (ref. 5) provides acceptable predictions of nonequilibrium and 
hypersonic boundary layers with cold walls and with blowing for a few investigations 
where experimental data were available. In a review of prediction methods for com- 
pressible turbulent boundary layers (ref. 6) the author concluded that finite -difference 
procedures offer the greatest potential for future uses. However, development of these 
methods and extension of their range of applicability requires detailed experimental 
measurements of both mean and fluctuating flow properties for a wide range of conditions. 

The lack of reliable and detailed data for compressible boundary layers has led 
to the development of transformation methods such as those presented in references 7 
and 8. These methods transform a given compressible -flow problem to a corresponding 
incompressible -flow problem, where the data are more accurate and plentiful. Exper- 
imental measurements of both pitot pressures and stagnation temperatures in compress- 
ible boundary layers are required to assess the accuracy and range of applicability of the 
transformations. 

Experimental data for hypersonic turbulent boundary layers are available in ref- 
erences 9 to 29. References 9 to 11 are for nearly adiabatic wall conditions. Data for 
conditions of large heat transfer with ratios of wall temperature to settling-chamber 
temperature T w /T 0 < 0.4 are reported in references 12 to 19. Detailed profile data 
(including measured total temperatures) for large heat -transfer conditions are available 
in references 12 to 16. Detailed profile data for 0.4 < T w /T 0 <0.5 are also available 
in references 21, 23, 24, and 25. 

Experimental studies of the boundary layer on the walls of a hypersonic nozzle are 
useful not only for direct applications to nozzle flow and nozzle design problems but also 
in the assessment and development of theoretical methods. In particular, these studies 
are of interest where upstream effects may be present (ref. 5). Also, since the boundary 
layers on nozzle walls are generally much thicker than on typical wind-tunnel models, 
more accurate and detailed measurements can be obtained in the wall boundary layer than 
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on models at comparable stream conditions. For Mach numbers greater than about 12, 
the boundary layers on wind-tunnel models will usually be laminar, whereas the nozzle - 
wall boundary layers will often be turbulent because of the much larger momentum thick- 
ness Reynolds numbers on the nozzle wall. 

The results reported herein are based on measurements of pitot and static pres- 
sure and total temperature in the turbulent boundary layer on the wall of a contoured noz - 
zle at a nominal stream Mach number of 19.5 with large heat -transfer rates to the wall 
(t w /T d ~0.2). The boundary layers were about 10.2 centimeters (4 in.) thick and the 
data were obtained with test duration times on the order of 2 hours. Wall skin friction 
was inferred from the slope of the Mach number profiles at the wall. Surface heat trans- 
fer was estimated from the slope of the total -temperature profile at the surface. No 
direct measurements of wall skin friction or heat transfer were made. 

Comparisons are made of the present test results with previous data at lower Mach 
numbers in air and nitrogen and with data in helium for Mach numbers up to 37. The 
present test data are also compared with results from several theoretical solutions by 
the finite -difference method of reference 5. By this means, the possible effects of lami- 
narization and turbulent Prandtl number have been investigated. Some preliminary 
results from the present investigation were compared with theoretical predictions in 
reference 5. For the final results presented herein, the pitot pressure data have been 
corrected for viscous and rarefaction effects, the corrections for heat conduction and 
radiation losses from the temperature probe have been improved, and the free -stream 
conditions have been modified slightly due to new measurements of static pressure. 

Some comparisons between the preliminary data as presented in reference 5 and the final 
corrected data are included in this report. Complete tabulations of the original uncor- 
rected data and the final corrected data are also given in the present report. 

SYMBOLS 

A constant in static-pressure distribution, equation (1) 

A* constant in equation (11), 26 except for wall blowing cases 

A,B,C functions of time in relation for wire temperature, equation (A2) 
a speed of sound 

D exponent in relation between N, q , and Tw/To 
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Cf 


d 

E 


specific heat of wire 


local skin friction coefficient, ■v/(i PeV 


specific heat at constant pressure, generally assumed constant except in 
determination of kj from Np r ^ 

diameter of wire 

voltage drop over specified length of wire 


voltage drop per unit length 


F 

f 

H 




h 

I 

K 

k 

L 

l 

M 

N 

N Nu,t 


velocity profile, u/u e 

mixing length function, equation (11) 

probe outside height 

incompressible form factor, C (1 - F)dy / f F(1 - F)dy 

Jo No 

probe inside height 


Me du e 


Pe u e 


2 dx 


electric current 
laminarization parameter, 
thermal conductivity 
exposed wire length 
mixing length 
Mach number 


reciprocal of exponent in power -law velocity relation, equation (15) 
Nusselt number, qd j (T aw - T w jk^ 
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N 


N Re,S 

N R e ,0 

% 

P 

q 

R 


T 

t 

u 


Prandtl number, Cpp/k 

turbulent Prandtl number based on total-temperature gradients 


Pr 

Npr,T 

Nr e/ m unit Reynolds number, pu/p, per meter 

Nr e/ft 
NRe,x 
N Re, 1 


unit Reynolds number, pu/p, per foot 


local Reynolds number based on streamwise distance from flow origin 


pitot -probe Reynolds number based on local conditions ahead of probe 
shock, PjUjH I Hi 

boundary -layer -thickness Reynolds number, p e u e 6/p e 
momentum -thickness Reynolds number, p e u e 0/p e 
Stanton number, 
pressure 


q w 


Pe u e c p(^aw " T w) 


heat -transfer rate 

ideal gas constant, 296.5 m^/sec^-K (l770 ft2/sec2-°R) used for present 
data 

radius from axis of symmetry 
recovery factor, T aw /Tj. 

absolute temperature 


time 


velocity parallel to wall in direction of free -stream flow 


velocity normal to wall 


W 


probe outside width 
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x axial distance from nozzle throat, or distance along thermocouple wire 

y normal distance from wall 

j3 local wall angle with respect to nozzle center line 

r intermittency function 

y ratio of specific heats, 1.4 used for present data 

6 boundary -layer thickness at point where M ~ 0.995M e unless otherwise 

noted 

6* displacement thickness, equation (5) 

e eddy viscosity or total hemispherical emissivity 

6 momentum thickness, equation (6) 

Q normalized total-temperature parameter, - T w j^ (T^ e - T w ^ 

p. viscosity 

p density 

a Stephan-Boltzmann constant 

t shear stress 

o /p w T p 1 

v viscous interaction parameter, M P — — 

x 6 W T w N R6jX 



Subscripts: 


A alumel 

a ambient 

av average 
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aw adiabatic temperature of wall or of bare -wire probe 

C chromel 

c longitudinal curvature 

e at "edge” of boundary layer 

f far-wall region 

F.M. free molecule 

i affected by wall interference 

M thermocouple junction 

max maximum 

meas measured 

n near-wall region 

o settling chamber 

p probe tip 

ref reference continuum value 

S wire temperature at x = ±L/2 

s viscous sublayer edge 

T turbulent or based on total-temperature profiles 

t local total or stagnation conditions 

w surface; or wire when applied to bare -wire probe 



X 


x -direction 


y y -direction 

1 local undisturbed flow (ahead of normal shock caused by probe) 

2 behind normal shock 

9 based on momentum thickness 

A primed symbol denotes a time fluctuating quantity and a superscript bar, except 
where otherwise noted, indicates a time mean value or correlation of the product of 
fluctuating quantities. 


APPARATUS AND TEST PROCEDURES 
Wind Tunnel 

The experimental data were obtained in the Langley hypersonic nitrogen tunnel. 

A detailed description of the facility and its operation is given in reference 30. This 
facility has a water-cooled axisymmetric contoured nozzle with an exit diameter of 
40.6 centimeters (16 in.) and an open-jet test section. The nozzle is about 3.2 meters 
(10.5 ft) in length (at the exit, x = 2.248 meters (88.5 in.)) and has wall static-pressure 
ports approximately every 0.3048 meter (1 ft) along the nozzle length. The facility can 
be operated for sustained periods of time (run times in excess of 2 hours) at stagnation 
pressures up to 6849 N/cm 2 (10 000 lb/in 2 ) and stagnation temperatures up to 2220 K 
(4000° R). The test gas is nitrogen with less than 5 parts per million of oxygen. 

The nominal free -stream conditions for the present tests were at a Mach number 
of 19.5 and over a range of Reynolds numbers per meter of 1.9 x 10® to 3.0 X 10® 

(5.8 x 10® to 9.2 x 10® per ft). The nominal total temperature was 1780 K (3200° R) 
and the wall -to -stagnation temperature ratio was approximately 0.17 at the boundary- 
layer-survey station which is 2.083 meters (82 in.) downstream of the throat. The var- 
iation in Mach number across the free -stream flow region is shown in reference 30. 

The axial variation in Mach number was found to be about 1.64 per meter (0.5 per ft) 
in the vicinity of the present survey station. 

Instrumentation 

The pressure in the settling chamber was measured with a strain-gage transducer. 
The nominal accuracy of this gage, according to the manufacturer, is ±1 percent of full 
scale which is 10.3 x 10 7 N/m 2 (15 000 psia). The lowest stagnation pressure used 
in the present investigation was 2.7 x 10 7 N/m 2 (4000 psia) and by careful calibration 
procedures, the accuracy for this pressure range was held to ±2 percent of the indi- 
cated reading. 
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The stagnation temperature was measured in the settling chamber with a platinum - 
platinum 13 percent rhodium thermocouple attached to the tip of a 0.32-centimeter 
(1/8 -in.) diameter, swaged, thin -wall tube that enclosed the thermocouple wire and high 
temperature insulation. The tube, or total -temperature probe, was mounted within the 
settling chamber just upstream of the approach section to the nozzle. (See ref. 30.) 
Corrections for heat conduction losses were neglected because of the construction and 
location of the probe. Estimated corrections due to radiation were from 3 to 5 percent 
of the measured probe temperature depending on the temperature assumed for surrounding 
surfaces. However, no radiation corrections were used because of uncertainties in the 
surrounding surface temperatures. 

An independent check of the stagnation temperature was obtained from a heat bal- 
ance based on the known mass-flow rate, the measured power input to the tungsten heater 
element, the measured air temperature at the inlet of the settling chamber, and estimated 
radiation and conduction losses. Values of total temperature obtained from this heat bal- 
ance were generally about 3 percent higher than the measured values. However, because 
of uncertainties in radiation and conduction losses, temperatures shown or used in the 
present report are the measured values with no corrections. This procedure was justified 
partly on the basis of the aforementioned uncertainties in radiation and conduction effects 
and partly due to the relatively small effect of a nominal increase in T 0 on the normal- 
ized total-temperature profiles. Thus, the maximum effect of a 5-percent increase in T 0 
on the values of Tt/T 0 was about 6 percent. 

Sketches of the boundary -layer probes and water-cooled strut support are given in 
figure 1. Most of the pitot -pressure data reported herein were obtained with the flattened 
stainless-steel tube mounted on the end of the strut next to the wall as indicated in the 
figure. This flattened tube had external dimensions at the mouth of 0.381 centimeter 
(0.15 in.) wide by 0.028 centimeter (0.011 in.) high. The total-temperature probe was a 
chromel-alumel bare -wire thermocouple of 0.025 -centimeter (0.01 -in.) diameter placed 
normal to the flow. The length -to -diameter ratio of the portion of the thermocouple nor- 
mal to the flow was 25. 

The location of the survey probes normal to the wall surface and in the streamwise 
direction was controlled by a precision boundary -layer -survey mechanism accurate to 
0.0025 centimeter (0.001 in.). The probes were traversed normal to the nozzle wall by 
remote control and a low voltage contact indicator stopped the movement at the wall sur- 
face, The data were continuously monitored during the tests and were not used until the 
pressure or temperature readings reached a constant level for a selected position in the 
boundary layer. 

Both the static and the pitot pressures were measured with absolute pressure trans- 
ducers. The basic pressure sensor is a high precision, stable capacitive bridge, the 
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Nozzle wall 



Side view of strut with pitot tubes 


, Figure l.- Schematic sketch of probe and strut. Dimensions in inches (parenthetically in centimeters). 



variable element of which is a thin, prestressed metallic diaphragm. The diaphragm is 
positioned between fixed capacitor plates, and its deflection is directly proportional to 
the magnitude of the pressure input. A regulated carrier voltage is applied across the 
capacitive bridge and the deflection of the diaphragm modulates the carrier voltage to 
provide an electrical output. The instruments used during this investigation have a range 
selection feature that provided seven ranges with full-scale readout on each range. Static 
pressures were measured with instruments which have a 0- to 50-torr scale on their 
highest range and a 0- to 0.5-torr scale on their lowest range. Pitot pressures were 
measured with instruments of 100-tor r full scale on the highest range and 1-torr scale 
on the lowest range. Continuing checks of the calibrations indicate that accuracies of 
±2 percent, or better, of full-scale readout were obtained. 

Test Conditions 

A summary of the nominal test conditions for the experimental program is provided 
in table I. The temperature surveys were obtained during different runs than the pitot 
surveys, and the corresponding test conditions were somewhat different as indicated in 
the table. For the duration of a run (up to about 2 hours), the nominal settling -chamber 
conditions varied a maximum of about 9 percent. Hence, all boundary -layer probe data 
have been normalized by the appropriate stagnation values recorded at the same time as 
the probe data. 

The nominal free -stream conditions at the survey station (x = 2.083 meters (82 in.)) 
where r w = 19.70 centimeters (7.758 in.) as given in table I were calculated by assuming 
isentropic expansion to the measured pitot pressures of an ideal gas with y = 1.4. The 
ideal gas properties of nitrogen (R = 296.5 m2/sec2-K (1770 ft2/sec2-°R)) have been used 
in all boundary -layer data-reduction procedures in the present investigation. The use of 
real gas properties, as obtained, for example, from reference 31, would affect computed 
flow properties by, at most, only a few percent. In view of other uncertainties in correc- 
tions for boundary -layer probe data, the use of real gas properties was not warranted for 
the present range of conditions. The viscosity for nitrogen was computed from equa- 
tion (44) of reference 32. 

TABLE I.- NOMINAL STAGNATION AND TEST CONDITIONS AT x = 2.083 METERS (82 IN.) AND T w = 300 K (540° R) 
(First and second sets of values for any given condition are for pitot-pressure and stagnation-temperature runs, respectively] 


Symbol 

Condition 

1 

0 

p< 

5 

p t,2,e/Po 

M e 

N Re/ m 

NRe/ft 

p e /p 0 

Pe/P o 

Pw/Pq 


K 

°R 

N/cm 2 

lb/in 2 

O 

1 

1680 

3020 

2970 

4300 

0.1291 X 10“ 3 

19.28 

2.02 x 10 6 

6.18 x 10 5 

0.269 x 10' 6 

0.203 X 10" 4 

0.391 x KT 6 



1695 

3050 

2830 

4100 



1.90 

5.78 




□ 

2 

1780 

3200 

4310 

6250 

.1245 

19.42 

2.62 

7.98 

.256 

.1958 

.389 



1695 

3050 

4380 

6350 



2.90 

8.84 




<0 

3 

1835 

3300 

5550 

8050 

.1175 

19.65 

2.98 

9.08 

.236 

.1848 

.390 



1835 

3300 

5590 

8100 



3.00 

9.16 
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PROBE ERRORS AND CORRECTION TECHNIQUES 


A complete tabulation of the uncorrected pitot pressures and total temperatures 
measured in the boundary layer and normalized with their respective settling-chamber 
values is given in table n. Also included in table II for comparison are the corresponding 
corrected values. The final corrected data are also given in table IH together with the 
resulting profiles of Mach number, velocity, static temperatures, and total temperature. 
The various sources of errors in probe data and the procedures used herein to determine 
the corrections will now be discussed in detail. 

TABLE H.- UNCORRECTED AND CORRECTED VALUES OF PITOT PRESSURES AND TOTAL TEMPERATURES 

AT x ~ 2.083 METERS (82 IN.) 


(a) Condition 1, M e = 19.28 




I 

> ~ 2970 N/cm 2 (4300 lb/in 2 ) 





P 0 S 

2830 N/cm 2 (4100 lb/in 2 ) 


y 

p t,2/ p o 



p t, 

2/ p o 



7 

Tt/T 0 

1 

"o 

cm 

in. 

Uncorrected 

Corrected 

cm 

in. 

Uncorrected 

Corrected 

cm 

in. 

Uncorrected 

Corrected 

K 

°R 

0.2743 

0.108 

1.059 x 10~ 6 

0.580 x 10-6 

10.180 

4.008 

125.2 x 10* 6 

126.37 x 10“6 

0.2019 

0.0795 

0.226 

0.225 

1678 

3020 

.7823 

.308 

1.261 

.639 

10.434 

4.108 

127.1 

127.10 

.7099 

.2795 

.227 

.349 

1689 

3040 

1.0363 

.408 

1.419 

.690 

10.688 

4.208 

128.3 

128.30 

1.2179 

.4795 

.309 

.423 

1694 

3050 

1.2903 

.508 

1.820 

.866 

10.942 

4.308 

128.5 

128.50 

1.7259 

.6795 

.347 

.515 

1706 

3070 

1.5443 

.608 

2.316 

1.092 

11.196 

4.408 

128.5 

129.10 

2.2339 

.8795 

.392 

.625 

1706 

3070 

1.7 983 

.708 

2.740 

1.282 

11.450 

4.508 

131.0 

131.0 


2.7419 

1.0795 

.4375 

.735 

1717 

3090 

2.0523 

.808 

3.927 

1.853 

11.704 

4.608 

131.0 

No correction 

3.2499 

1.2795 

.489 

.755 

1667 

3000 

2.3063 

.908 

4.996 

2.393 

11.958 

4.708 

131.0 



3.7579 

1.4795 

.521 

.772 

1700 

3060 

2.5603 

1.008 

6.159 

2.994 

12.212 

4.808 

132.0 



4.2695 

1.6795 

.549 

.795 

1694 

3050 

2.8143 

1.108 

7.368 

3.644 

12.466 

4.908 

131.0 



4.7739 

1.8795 

.574 

.792 

1689 

3040 

3.0683 

1.208 

9.016 

4.566 

12.720 

5.008 

131.0 



5.2819 

2.0795 

.591 

.805 

1689 

3040 

3.5763 

1.408 

13.380 

7.378 

12.974 

5.108 

130.0 



5.7899 

2.2795 

.608 

.825 

1706 

3070 

3.8303 

1.508 

16.370 

9.466 

13.228 

5.208 

129.9 



6.2979 

2.4795 

.623 

.862 

1711 

3080 

4.0843 

1.608 

19.920 

11.918 

13.736 

5.408 

124.5 



6.8059 

2.6795 

.636 

.885 

1722 

3100 

4.3383 

1.708 

23.250 

14.143 

13.990 

5.508 

124.5 



7.3139 

2.8795 

.649 

.899 

1739 

3130 

4.5923 

1.808 

25.280 

15.347 

14.244 

5.608 

124.7 



7.8219 

3.0795 

.653 

.919 

1761 

3170 

4.8463 

1.908 

29.68 

18.008 

14.498 

5.708 

i 127.0 



8.3299 

3.2795 

.661 

.932 

1772 

3190 

5.1003 

2.008 

33.28 

21.965 

14.752 

5.808 

128.5 



8.8379 

3.4795 

.672 

.950 

1778 

3200 

5.3543 

2.108 

37.74 

27.570 

15.006 

5.908 

129.0 



9.3459 

3.6795 

.677 

.959 

1789 

3220 

5.6083 

2.208 

41.49 

32.803 

15.260 

6.008 

126.5 



9.8539 

3.8795 

.684 

.971 

1794 

3230 

5.8623 

2.308 

45.51 

38.315 

15.551 

6.108 

128.5 



10.362 

4.0795 

.690 

.980 

1794 

3230 

6.1163 

2.408 

49.84 

43.734 

15.768 

6.208 

130.0 



10.87 

4.2795 

.694 

.990 

1794 

3230 

6.3703 

2.508 

54.34 

49.124 

16.022 

6.308 

132.0 



11.378 

4.4795 

.698 

1.00 

1794 

3230 

6.6243 

2.608 

58.31 

53.796 

16.276 

6.40J 

133.0 



11.863 

4.6705 

.700 

1.00 

1794 

3230 

6.8783 

2.708 

63.40 

59.598 

16.530 

6.508 

136.5 



12.394 

4.8795 

.700 

1.00 

1794 

3230 

7.3863 

2.908 

74.58 

72.298 

16.784 

6.608 

138.0 









7.6403 

3.008 

80.4 

78.648 

17.038 

6.708 

139.0 









7.8943 

3.108 

85.76 

84.487 

17,546 

6.908 

138.5 









8.1483 

3.208 

90.96 

90.115 

17.800 

7.008 

136.0 

\ 

t 







8.4023 

3.308 

95.23 

94.741 












8.6563 

1 3.408 

99.8 

99.599 












8.9103 

3.508 

104.0 

104.11 












9.1643 

3.608 

109.6 

110.15 












9.4183 

3.708 

114.3 

115.06 












9.6723 

3.808 

119.6 

120.57 












9.9263 

3.908 

123.1 

124.29 
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Probe Interference Effects 


The question always arises as to whether or not the flow field is affected by the 
presence of survey probes. Various investigators have recognized the wall- interference 
problem and considered possible effects on pitot- or temperature -probe measurements 
(see, for example, refs. 19 and 33 to 35). Most investigators eliminate or discount data 

TABLE n.- UNCORRECTED AND CORRECTED VALUES OF PITOT PRESSURES AND TOTAL TEMPERATURES 

AT x = 2.083 METERS (82 IN.) - Continued 

(b) Condition 2, M e = 19.42 


y 

p Q = 4310 N/cm 2 (6250 lb/in 2 ) 



Po = 

4380 N/cm 2 (6350 lb/in 2 ) 




p t,2/Po 


y 


Tt/T 0 

T 

o 

cm 

in. 

Uncorrected 

Corrected 

cm 

in. 

Uncorrected 

Corrected 

K 

OR 

0.2743 

0.108 

0.765 X 10- 6 

0.4692 x 10“6 

0.2019 

0.0795 

0.226 

0.253 

1700 

3060 

.5283 

.208 

.9847 

.5462 

.6337 

.2495 

.274 

.325 

1672 

3010 

.7823 

.308 

1.354 

.6748 

1.1417 

.4495 

.335 

.438 

1689 

3040 

1.0363 

.408 

1.897 

.9282 

1.6497 

.6495 

.405 

.593 

1656 

2980 

1.2903 

.508 

2.439 

1.1834 

2.1577 

.8495 

.455 

.689 

1672 

3010 

1.5443 

.608 

3.412 

1.6651 

2.6657 

1.0495 

.520 

.751 

1644 

2960 

1.7983 

.708 

4.578 

2.2888 

3.1737 

1.2495 

.555 

.782 

1667 

3000 

2.0523 

.808 

5.903 

3.0163 

3.6817 

1.4495 

.585 

.787 

1667 

3000 

2.3063 

.908 

6.815 

3.5255 

4.1897 

1.6495 

.605 

.795 

1678 

3020 

2.5603 

1.008 

8.749 

4.7361 

4.6977 

1.8495 

.621 

.814 

1689 

3040 

2.8143 

1.108 

10.71 

6.0718 

5.2057 

2.0495 

.640 

.819 

1656 

2980 

3.0683 

1.208 

13.11 

7.8149 

5.7137 

2.2495 

.659 

.841 

1644 

2960 

3.3223 

1.308 

16.39 

10.306 


6.2217 

2.4495 

.680 

.882 

1694 

3050 

3.5763 

1.408 

20.05 

13.079 


6.7297 

2.6495 

.693 

.903 

1700 

3060 

3.8303 

1.508 

22.96 

15.142 


7.2377 

2.8495 

.710 

.931 

1744 

3140 

4.0843 

1.608 

26.52 

17.617 


7.7457 

3.0495 

.721 

.954 

1744 

3140 

4.3383 

1.708 

29.03 

20.658 


8.2537 

3.2495 

.730 

.966 

1755 

3160 

4.5923 

1.808 

32.80 

25.766 


8.7617 

3.4495 

.732 

.979 

1800 

3240 

4.8463 

1.908 

36.84 

31.378 


9.2697 

3.6495 

.743 

.999 

1767 

3180 

5.1003 

2.008 

42.03 

37.785 


10.794 

4.2495 

.750 

1.00 

1756 

3160 

5.3543 

2.108 

46.35 

42.999 


11.81 

4.6495 

.732 

1.00 

1783 

3210 

5.6083 

2.208 

50.07 

47.264 








5.8623 

2.308 

55.42 

53.411 








6.1163 

2.408 

59.50 

57.973 








6.3703 

2.508 

65.95 

65.017 








6.6243 

2.608 

70.86 

70.220 








6.8783 

2.708 

76.14 

75.891 








7.1323 

2.808 

81.49 

81.668 








7.3863 

2.908 

87.4i 

87.753 








7.6403 

3.008 

92.02 

92.556 








7.8943 

3.108 

96.99 

97.686 








8.1483 

3.208 

104.4 

105.34 








8.4023 

3.308 

109.8 

110.82 








8.6563 

3.408 

117.1 

118.41 








8.9103 

3.508 

119.5 

121.31 








9.1643 

3.608 

121.0 

122.93 








9.4183 

3.708 

123.2 

124.50 








9.6723 

3.808 

124.0 

No correction 







9.9263 

3.908 

124.2 









10.18 

4.008 

124.6 









10.434 

4.108 

126.2 









10.942 

4.308 

124.5 









11.45 

4.508 

122.0 

> 
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obtained in the region of wall interference and extrapolate trends indicated by remaining 
data to the wall. 


In order to determine possible flow -field disturbances caused by the probes in the 
present tests, the wall pressure was monitored continuously as the survey probes were 
traversed away from the nozzle wall. Typical results are shown in figure 2 for surveys 
with the total -temperature probe and with the pitot -pres sure probe at x = 2.083 meters 
(82 in.). The probe tips were located opposite the downstream edge of the nozzle-wall 
static -pressure orifice as indicated in the figure. Measured values of p w ^/p Q are 

TABLE II.- UNCORRECTED AND CORRECTED VALUES OF PITOT PRESSURES AND TOTAL TEMPERATURES 

AT x= 2.083 METERS {82 IN.) - Concluded 

(c) Condition 3, M e = 19.65 




p Q = 5550 N/cm 

2 (8050 lb/in 2 ) 




Pq * 

5590 N/cm 2 (8100 lb/in 2 ) 


y 

P t,2/Po 


5 


Px « /p 

y 

Tt/' 

r 0 


Co 

cm 

in. 

Uncorrected 

Corrected 

cm 

in. 

2/ *o 

cm 

in. 

Uncorrected 

Corrected 

K 

°R 

0.2743 

0.108 

0.735 X 10" 6 

0.4591 x 10~6 

10.434 

4.108 

127.0 x 10-6 

0.2019 

0.0795 

0.208 

0.227 

1833 

3300 

.5283 

.208 

1.12 

.5992 

10.688 

4.208 

125.5 

.7099 

.2795 

.290 

.385 

1833 

3300 

.7823 

.308 

1.68 

.8286 

10.942 

4.308 

122.0 

.9639 

.3795 

.361 

.501 

1822 

3280 

1.0363 

.408 

2,44 

1.2026 

11.196 

4.408 

119.0 

1.4719 

.5795 

.428 

.641 

1850 

3330 

1.2903 

.508 

3,41 

1.6811 

11.45 

4.508 

118.0 

1.9799 

.7795 

.500 

.762 

1844 

3320 

1.5443 

.608 

4.84 

2.4515 

11.704 

4.608 

113.0 

2.4879 

.9795 

.541 

.786 

1850 

3330 

1.7983 

.708 

7.11 

3.8144 

11,958 

4.708 

112.0 

2.9959 

1.1795 

.572 

.795 

1850 

3330 

2.0523 

.808 

9.14 

5.2117 

12.212 

4.808 

111.0 

3.5039 

1.3795 

.600 

.808 

1844 

3320 

2.3063 

.908 

11.85 

7.1397 

12.466 

4.908 

109.0 

4.0119 

1.5795 

.621 

.823 

1844 

3320 

2.5603 

1.008 

14.63 

9.4464 

12.72 

5.008 

107.5 

4.5199 

1.7795 

.636 

.834 

1844 

3320 

2.8143 

1.108 

17.85 

11.946 

12.974 

5,108 

104,5 

5.0279 

1.9795 

.656 

.858 

1844 

3320 

3.0682 

1.208 

21.0 

14.396 

13.228 

5.208 

103.5 

5.5359 

2.1795 

.673 

oo 

00 

1844 

3320 

3.3223 

1.308 

24.65 

16.927 

13.482 

5.308 

103.0 

6.0439 

2.3795 

.686 

.903 

1844 

3320 

3.5763 

1.408 

28.35 

20.962 

13.736 

5.408 

105.0 

6.6519 

2.5795 

.707 

.929 

1844 

3320 

3.8303 

1.508 

32.35 

26.422 

13.99 

5.508 

105.0 

7.0599 

2.7795 

.721 

.952 

1844 

3320 

4.0843 

1.608 

35.45 

32.161 

14.244 

5.608 

103.5 

7.5679 

2.9795 

.736 

.974 

1844 

3320 

4.3383 

1.708 

41.2 

37.336 

14.498 

5.708 

101.5 

8.0759 

3.1795 

.744 

.984 

1844 

3320 

4.5923 

1.808 

46.1 

43.167 

14.752 

5.808 

102.0 

8.5839 

3.3795 

.751 

.999 

1844 

3320 

4.8463 

1.908 

51.2 

49.016 

15.006 

5.908 

103.5 

9.0919 

3.5795 

.759 

1.000 

1844 

3320 

5.1003 

2.008 

55.7 

54.360 

15.26 

6.008 

105.0 

9.5999 

3.7795 

.756 

1.000 

1844 

3320 

5.3543 

2.108 

61.1 

60.346 

15.514 

6.108 

101.5 

10.108 

3.9795 

.761 

1.000 

1844 

3320 

5.6083 

| 2.208 

! 66.0 

65.786 

15.768 

6.208 

99.5 

10.616 

4.1795 

.759 

1.000 

1833 

3300 

5.8623 

2.308 

71.4 

71.497 

16.022 

6.308 

99.5 

11.124 

4.3795 

.756 

1.000 

1833 

3300 

6.1163 

2.408 

77.0 

76.841 

16.276 

6.408 

99.5 

11.632 

4.5795 

.753 

1.000 

1839 

3310 

6.3703 

2.508 

82.9 

83.585 

16.53 

6.508 

100.0 

12.648 

4.9795 

.746 

1.000 

1839 

3310 

6.6243 

2.608 

88.2 

89.125 

16.784 

6.608 

99.9 







6.8783 

2.708 

94.4 

95.385 

17.038 

6.708 

95.5 







7.1323 

2.808 

100.0 

101.03 


17.292 

6.808 

94.5 







7.3863 

2.908 

104.0 

105.43 


17.546 

6.908 

93.0 







7.6493 

3.008 

108.0 

109.44 











7.8943 

3.108 

112.8 

114.55 











8.1483 

3.208 

116.7 

118.42 











8.4023 

3.308 

119.2 

121.49 











8.6563 

3.408 

119.7 

119.50 











8.9103 

3.508 

121.1 

No correction 










9.1643 

3.608 

123.9 












.4183 

3.708 

125.8 












9.6723 

3.808 

126.0 












9.9263 

3.908 

127.5 












10.18 

4.008 

127.5 


f 
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TABLE m.- EXPERIMENTAL BOUNDARY- LAYER PROFILES WITH 
ALL CORRECTIONS INCLUDED 


(a) Condition 1: 6 = 11.2 cm (4.4 in.); M e = 19.28; 6 = 0.178 cm (0.070 in.); 

6* = 5.92 cm (2.33 in.); A = 1.0 



M[M e 


0.0373 

.0446 

.0581 

.0804 

.1089 

.1426 

.1805 

.2256 

.2723 

.3294 

.4072 

.4876 

.5576 

.6276 

.6950 

.7573 

.8117 

.8640 

.9150 

.9611 

.9865 

.9984 

1.000 


T/T e 

15.10 
21.15 

24.10 
24.72 
23.78 
21.03 
16.98 
12.97 

9.74 

6.95 

4.67 

3.34 

2.65 

2.16 

1.81 

1.55 

1.39 

1.25 

1.14 

1.05 

1.01 

.99 

1.00 


T t /T 0 

0.225 

.349 

.423 

.515 

.625 

.735 

.755 

.772 

.795 

.792 

.805 

.825 

.862 

.885 

.899 

.919 

.932 

.950 

.959 

.971 

.980 

.990 

1.00 


u/u, 


e 

0.1451 

.2063 

.2847 

.3993 

.5306 

.6535 

.7429 

.8116 

.8489 

.8676 

.8786 

.8907 

.9071 

.9210 

.9348 

.9425 

.9548 

.9658 

.9740 

.9822 

.9871 

.9932 

1.000 


0.072 

.206 

.300 

.408 

.547 

.679 

.763 

.754 

.768 

.773 

.791 

.835 

.864 

.872 

.895 

.925 

.940 

.953 

.964 

.984 

.999 

1.00 


plotted against the distance of the probe from the wall. Wall static pressure was 
increased by a maximum of about 17 percent of the interference -free value p w (value 
at which the measured wall static pressure remains nearly constant with increasing dis- 
tance of probe from wall) when the pitot probe was near the wall. Estimates of maximum 
interference effects when the local flow is supersonic were obtained by assuming that 
oblique shocks form in front of the pitot tube (see refs. 34 and 35) and that the static - 
pressure ratio across these shocks is given by p w ^ /p w as obtained from figure 2. 
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TABLE m.- EXPERIMENTAL BOUNDARY-LAYER PROFILES WITH 
ALL CORRECTIONS INCLUDED - Continued 


(b) Condition 2: 6 = 9.4 cm (3.7 in.); M e = 19.42; 6 = 0.173 cm (0.068 in.); 

5* = 4.62 cm (1.82 in.); A = 1.0 


y/6 

P/Pt,2 

M/M e 

T/T e 

H 

c+ 

i-3 

o 

CD 

9 

0.0213 

859.6 X 10 -3 

0.0242 

17.73 

0.253 

0.1023 

0.0932 

.0671 

643.0 

.0422 

20.62 

.325 

.1924 

.1775 

.1211 

359.9 

.0680 

23.60 

.438 

.3311 

.316 

.1750 

187.3 

.0999 

24.09 

.593 

.4921 

.502 

.2290 

103.6 

.1375 

20.42 

.689 

.6236 

.621 

.2828 

47.47 

.2060 

12.63 

.751 

.7347 

.695 

.3368 

29.45 

.2626 

9.03 

.782 

.7923 

.735 

.3906 

18.40 

.3332 

6.01 

.787 

.8200 

.740 

.4445 

12.22 

.4094 

4.20 

.795 

.8424 

.751 

.4983 

8.57 

.4892 

3.10 

.814 

.8649 

.774 

.5520 

6.34 

.5690 

2.29 

.819 

.8650 

.796 

.6062 

5.00 

.6411 

1.86 

.841 

.8767 

.804 

.6600 

4.13 

.7055 

1.67 

.882 

.9142 

.856 

.7140 

3.53 

.7631 

1.47 

.903 

.9281 

.884 

.7680 

3.05 

.8213 

1.35 

.931 

.9561 

.918 

.8220 

2.66 

.8795 

1.21 

.954 

.9689 

.944 

.8740 

2.37 

.9310 

1.09 

.966 

.9789 

.958 

.9300 

2.17 

.9743 

1.04 

.979 

.998 

.975 

.9840 

2.72 

.9964 

1.01 

.999 

.999 

.998 

1.000 

2.09 

1.000 

1.00 

1.00 

1.00 

1.00 


The Mach number distribution downstream of the oblique shocks was obtained from the 
Rayleigh pitot -pressure formula using the measured pitot pressure (corrected for vis- 
cous and rarefaction effects) and the disturbed wall static pressure p w j, also from fig- 
ure 2. The Mach number distribution upstream of the assumed oblique shocks could then 
be computed and these values were, at most, 0.07 larger than the downstream values. 
Hence, no corrections for wall-interference effects have been applied to measured values 
of pitot pressure and stagnation temper attire, and the interference -free values of wall 
static pressures p w have been used in all data-reduction procedures which are dis- 
cussed in subsequent sections herein. 
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TABLE m.- EXPERIMENTAL BOUNDARY-LAYER PROFILES WITH 
ALL CORRECTIONS INCLUDED - Concluded 

(c) Condition 3: 5 = 8.65 cm (3.4 in.); M e = 19.65; 0 = 0.155 cm (0.061 in.); 

5 * =4.31 cm (1.70 in.); A = 1.0 


y/6 

P/Pt,2 

M/M e 

T/T e 

T t/T 0 

u/u e 

6 

0.0234 

859.6 x 10~ 3 

0.0239 

16.94 

0.227 

0.0982 

0.0756 

.0822 

498.0 

.0534 

25.30 

.385 

.2467 

.205 

.1116 

343.5 

.0692 

28.36 

.501 

.3677 

.403 

.1704 

160.8 

.1074 

26.66 

.641 

.5531 

.571 

.2292 

71.92 

.1643 

19.38 

.762 

.7218 

.716 

.2881 

38.52 

.2264 

12.47 

.786 

.7977 

.744 

.3469 

23.64 

.2901 

8.35 

.795 

.8359 

.756 

.4057 

15.30 

.3613 

5.72 

.808 

.8619 

.773 

.4645 

10.45 

.4376 

4.09 

.823 

.8826 

.788 

.5233 

7.44 

.5191 

3.00 

.834 

.8969 

.802 

.5822 

5.71 

.5929 

2.39 

.858 

.9148 

.831 

.641 

4.55 

.6641 

1.97 

.881 

.9299 

.859 

.6998 

3.72 

.7328 

1.67 

.903 

.9441 

.883 

.7587 

3.15 

.7990 

1.45 

.929 

.9593 

.917 

.8175 

2.72 

.8600 

1.29 

.952 

.9725 

.944 

.8763 

2.40 

.9160 

1.16 

.974 

.9848 

.970 

.9351 

2.15 

.9669 

1.06 

.984 

.9907 

.981 

.9939 

2.04 

.995 

1.02 

.999 

.998 

.998 

1.000 

2.05 

1.00 

1.00 

1.000 

1.000 

1.00 


When the local flow is subsonic and p w is used with measured pitot pressures to 
compute Mach number, the only wall -interference effect would be caused by displacement 
of streamlines away from the wall due to local flow separation as the probe approached 
the wall. The subsonic data should then be plotted at smaller values of y than indicated. 
Since the magnitude of this effect is probably small for the present conditions, no correc- 
tions have been applied. 


Pitot-Pressure Corrections 

Effect of thermal diffusion . - The pressure-survey data were obtained for nominal 
stagnation temperatures of 1670 to 1830 K (3000° to 3300° R). Because of the resulting 
high temperatures near the opening of the pitot tubes, the pressure measurements could 
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be affected by thermal diffusion since the probe strut was water cooled. An uncooled 
length of pitot tube extended about 3.2 centimeters (1.25 in.) ahead of the strut support. 

In order to compute corrections for this effect, the temperature distribution along the 
pitot tube would have to be known. Even if this temperature distribution were known, 
accurate corrections could not be computed for the present tests because of the complex 
geometry of the tip of the pitot tube. (See fig. 1.) Estimates of the thermal -diffusion 
effect were based on the assumptions that the temperature of the pitot -probe lip is the 
same as the measured total temperature (uncorrected) and that a mean-effective hydraulic 
radius for the probe tip region of 0.043 centimeter (0.017 in.) can be used in place of the 
radius of a circular tube. From the method of reference 36, the maximum error was 
then obtained as 4 percent when the temperature of the cold end of the tube was taken as 
300 K (540° R). Since these assumptions are probably conservative, no corrections for 
this effect have been used. 


V iscous and rarefied flow effects. - Studies of viscous and rarefied flow effects on 
impact pressure measurements have been reported in references 33 and 37 through 43. 

It was concluded (ref. 43) that the interpretation of impact pressure measurements at low 
density or low probe Reynolds numbers (as would be the case in the near -wall region) may 
be uncertain. Viscous effects on flattened probes were obtained over a limited Reynolds 
number range during the present investigation. Impact pressures were measured in the 
free stream with three tubes of different sizes and width to height ratios. The reference 
pressure was measured with a large pitot probe of 0.84-centimeter (0.33-in.) height. 

The results are compared with those of references 33 and 37 to 42 in figure 3. The sym- 
bol key used for the data in this figure and information on probe geometry are given in 
table IV. The faired lines represent the mean of data from all sources (short dashed 
lines are extrapolations) and have been used for the corrections applied herein. Alternate 
ways of plotting these data are illustrated in references 33 and 37 to 42; however, the 
present plot shows the Mach number effect explicitly and appeared to correlate the data 
better than Knudsen number or normal-shock Reynolds number PgUgHyV^.. 


Very little data are available for N Re j < 10 except at Mj < 2.0. Hence, for 
Nr 6 , i < 10, the faired lines used for pitot-pressure corrections are based on trends of 
available data and consideration of the appropriate free-molecule limiting values computed 
(formula (7-6), p. 702, ref. 44) with the assumption that the reflection coefficients are 
completely diffuse. The results of this calculation for the free-molecule limiting values 
are plotted on the insert in figure 3. The maximum effect occurs at Mach number 1.0 
and depends on the ratio of the probe temperature to stagnation temperature. The short- 
dashed-line curves in figure 3 represent arbitrary extrapolations to the free-molecule 
limits for Tp^T 0 ~ 1.0. 
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TABLE IV.- SYMBOL KEY FOR VISCOUS AND RAREFACTION CORRECTIONS TO PITOT PRESSURES AS SHOWN IN FIGURE 3 


Symbol 


M e 

H, 

cm 

h/H 

W/H 


(in.) 




Supersonic and hypersonic data 


o 

+ 

A 


1.7 to 2.0 


2.3 to 2.7 
2.8 to 3.4 


Reference tube 


1.23 to 1.99 


0.254 to 2.54 
(0.1 to 1.0) 

1 

1.52 to 2.54 
(0.6 to 1.0) 

0.254 to 2.54 
(0.1 to 1.0) 


4.4 to 4.6 

5.5 to 5.7 
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An iterative procedure (see appendixes A and B) was required to determine the 
local probe Reynolds number (Nj^ j = pjUjH/fij) in the present data-reduction procedure. 
Large corrections were obtained for pitot-pressure data in the near -wall region, and sig- 
nificant corrections were required up to 5.1 centimeters (2 in.) from the wall. Hence 
all data presented in table in have been corrected where required. In general, the cor- 
rections reduced the measured pitot pressure in the wall region since the rarefied flow 
effect caused the indicated measurements to read high in this region where N Re j < 100 
and M e i < 6.0. The magnitude of the corrections used through the boundary layer can 
be obtained for all runs from comparisons of values listed in table n for the uncorrected 
and corrected pitot pressures. Also, for the lowest pressure run (condition 1 in table I) 
where the corrections would be the largest, the final values of (p t 2)meas/( p t 2)ref are 
plotted against the corresponding values of Njj e ^ in figure 3 as the X symbols. For 
this condition the maximum correction ratio for the pitot pressure was 2.14 where the 
local probe Reynolds number was 0.45. The smallest probe Reynolds number for all 
three runs was 0.39 at y = 0.274 centimeter (0.108 in.). 

Effects of turbulent fluctuations . - The turbulent fluctuations in velocity and density 
will generally increase the local mean pitot pressures according to the analysis of ref- 
erence 45. The maximum increase for the present conditions was estimated to be about 
10 percent based on the data correlations and discussion of reference 45. This maximum 
increase in measured pitot pressure would occur in the vicinity of the viscous -sublayer 
edge which, in the present tests, was located about 1.91 centimeters (0.75 in.) from the 
wall. The corresponding maximum increase in Mach number would be about 6 percent. 

No attempt has been made to correct the pitot pressures for these effects because of 
uncertainties in the magnitude and distribution of the fluctuating quantities and of the 
second-order correlation coefficient (see ref. 45) for the present conditions. 

Another effect of turbulent fluctuations is manifested through the acoustic response 
of the measurement system (the probe, tubing, and transducer volume). This effect has 
been treated in reference 46 where it was shown that the mean pressures may be further 
increased due to the acoustic response of the system by amounts that are comparable to 
the increases noted in the previous paragraph. Since no general analysis of this effect 
is possible for supersonic turbulent flow, such effects are not considered further herein. 

Corrections for Radiation and Conduction Losses to Temperature Probe 

The data from the bare -wire chromel-alumel thermocouple used in the surveys 
were corrected for heat losses due to radiation and conduction. A detailed discussion of 
the method used and auxiliary data required for correcting the temperature survey data 
are given in appendixes A and B. Corrections due to conduction and radiation varied 
through the boundary layer and were 20 to 30 percent in the free stream. An indication 
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of the magnitude of the correction applied through the boundary layer may be obtained by 
comparing values from table II for the uncorrected and corrected total temperatures. 
Consideration of all sources of error in the data-reduction procedure for total tempera- 
ture (excluding any temperature -probe interference effects and any effects of fluctuating 
values) indicates that the final dimensionless profiles are correct to within approxi- 
mately ±10 percent. 


EXPERIMENTAL RESULTS 

Pressure and Temperature Measurements 

Pitot -pres sure profiles. - Figure 4 shows the variation in the ratio of pitot -to -total 
pressure with distance from the nozzle wall for the three test runs (table I). Corrections 
for viscous and rarefaction effects have been applied to the pitot data. The pitot -pressure 
data have been normalized by the settling-chamber pressure recorded at the same time 
to account for any slight variation in total pressure during the survey runs. Pitot- 
pressure data were obtained in the inviscid flow out to 17.8 centimeters (7 in.) from the 
nozzle wall for the test runs at the upper and lower total pressures. The degree of flow 
uniformity in the inviscid stream at the measuring station (x = 2.083 meters (82 in.)) is 
apparent from these data for y > 10.2 centimeters (4 in.). There is very little scatter 
in the experimental data, however, there is a slight peak in the pitot profiles as they 
join the inviscid core flow (8.6 cm < y < 12.7 cm (3.4 in. < y < 5.0 in.)). The nominal 
boundary-layer edge 6 has been chosen somewhat closer to the wall than the location 
of the peak pitot pressure for each test. While the resulting values of 6 cannot be 
specified accurately, they provide an approximate indication of boundary -layer thickness 
and in terms of y/5, the profiles tend to collapse to about the same distributions as will 
be shown in subsequent sections of this report. 

Previous investigations (see ref. 47) have shown that nozzle -wall turbulent boundary 
layers may differ from flat-plate boundary layers at similar conditions. These differ- 
ences would be due to upstream wall temperature and pressure gradients along the wall 
as well as normal pressure gradients, none of which would exist for uniform flat-plate 
flows. Any upstream wave system (resulting from an overexpansion at the throat, for 
example) would have some effect on the turbulent boundary -layer profiles. Also, such an 
upstream wave system could have some influence on the peak pitot pressure measured in 
the vicinity of the boundary -layer edge as seen in figure 4 at y * 10.2 centimeters (4 in.). 

In order to determine the effects of such wave systems or disturbances on the invis- 
cid flow field in the nozzle, a center-line pitot-pressure survey was made from about 
5.1 centimeters (2 in.) downstream of the nozzle throat to the exit (x = 2.248 meters 
(88.5 in.)). These data were used to compute the Mach number (with isentropic flow 
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Figure k.- Pitot -pres sure profiles corrected for viscous and rarefied flow effects 

(Uncorrected values given in table II.) 





assumed from the settling chamber) and the results are shown in figure 5. It is seen 
that disturbances exist for distances greater than about 1/3 of the nozzle length with a 
large overexpansion from x K 64 to 98 centimeters (24 to 37 in.). These disturbances 
tend to diminish further downstream as the x = 2. 083 -meter (82 -in.) station is 
approached. 

To investigate the origin of the overexpansion and subsequent disturbances in the 
nozzle, the inviscid flow field based on this center -line Mach number distribution was 
calculated by the method of characteristics as applied to nozzle flows in reference 48. 
Equilibrium real gas properties of nitrogen were used in this calculation with stagnation 
conditions taken as p Q = 4830 N/cm^ (7000 lb/in^) and T 0 = 1705 K (3070° R). (The 
data of fig. 5 were obtained with p Q = 3450 N/cm^ (5000 lb/in^); however, the computed 
results for this value of p Q would be nearly identical to those with p Q = 4830 N/cm^ 
(7000 lb/in2).) The results of this inviscid-flow-field calculation showed that the center - 
line overexpansion originates in a region on the wall just slightly downstream of the 
nozzle -wall inflection point. (This finding is consistent with the results of ref. 49, which 
indicates the overexpansion may be caused by the abrupt change in slope of the center - 
line Mach number distribution used in the original design calculation. For the present 
nozzle, this change in slope was at x* 92.60 centimeters (36.5 in.).) The lateral pitot- 
pressure distribution in the inviscid flow at x = 2.032 meters (80 in.) obtained from this 
calculation is shown in figure 4. The resulting wavy distribution is similar to the data, 
and indicates that the observed variations in the inviscid flow are caused by the overexpan- 
sion in the upstream core flow (fig. 5). These relatively small variations in the inviscid 
flow (fig. 4) are believed to have little effect on the wall boundary -layer development. 

Wall and free-stream static pressures .- Measured values of wall static -to -total 
pressure ratios and calculated values of free-stream static -to -total pressure ratios at 
x = 2.083 meters (82 in.) are shown in figure 6. The values of the calculated free-stream 
static pressures were obtained from the isentropic relation: 
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where M e is obtained from the isentropic relation between free-stream pitot pressure 
and settling -chamber pressure: 
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There is a significant difference between the measured wall pressures and calcu- 
lated static pressures shown in figure 6. These results are different from previous 
lower Mach number results (for example, refs. 12 and 19), where the agreement between 
calculated free-stream and measured wall static pressures was found to be good. How- 
ever, investigation of nozzle-wall boundary layers in helium at Mach numbers from 20 
to 47 (refs. 11 and 16) indicated that wall static pressures were larger than stream static 
pressures. 

Nozzle-wall pressures measured at several locations ahead of the survey station 
are also shown in figure 6. Included for comparison in the figure are ratios of stream 
static pressure to total pressure, shown as the solid line, obtained from the same inviscid- 
flow-field solution described in the preceding section. These values were obtained by 
interpolation in the computed characteristics mesh at the distances from the nozzle cen- 
ter line of (r w - 6) where 6 was the boundary -layer thickness computed by the finite- 
difference method of reference 5. (Solution no. 2 was used for this purpose. The com- 
puting procedure and inputs used for this solution are discussed in subsequent sections 
herein.) These values of static pressure tend to approach the measured wall pressures 
upstream of x = 102 centimeters (40 in.). The value at x = 2.083 meters (82 in.) is 
in good agreement with values obtained from measured pitot pressures at this station. 

The magnitude of the streamwise pressure gradient upstream of the survey station 
from x = 1.53 to 2.04 meters (60 to 80 in.) is approximately (3.54 X 10 -7)p o per meter 
or (o .9 X 10“8)p o per inch. The pressure drop across the boundary layer as indicated 
by the difference between the measured wall static values at x = 2.083 meters (82 in.) 
and the computed free-stream values (figs. 6 and 7) is about (0.15 x 10 - 6)p . The mag- 
nitude of the normal pressure gradient across the 10. 2 -centimeter (4 -in.) boundary layer 
is then roughly (l.57 x 10“®)p o per meter or (4 x 10-8)p o per inch or about four times 
larger than the streamwise pressure gradient. Hence, it might be expected that the nor- 
mal pressure gradient could have at least as large an effect on the boundary -layer struc- 
ture as the streamwise pressure gradient. The difference between the measured wall 
static pressures and computed free-stream static pressures may be caused by upstream 
disturbances, by turbulent velocity fluctuations, or by the inviscid-flow expansion that is 
present in all nozzle flows. If the latter mechanism affects the normal static -pressure 
gradients within the boundary layer, the pressure at the wall should be above the pressure 
near the edge of the boundary layer, as indicated by the results shown in figure 6. 

Static -pressure profi les.- Static -pressure distributions across the boundary layer 
obtained from measurements with a static -pressure probe and from various computations 
are shown in figure 7(a). Consider first the static -pres sure distribution (shown as the 
dashed line in fig. 7(a)) from the inviscid-flow -field calculation described previously. 

The aforementioned trend of increasing static pressure as the wall is approached is 
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Figure 7-- Concluded. 

indicated by the inviscid-flow calculation which was based on the center -line Mach num- 
ber distribution of figure 5. 

Also shown in figure 7(a) as an indication of wall -interference effects caused by the 
static probe are the values of the wall static pressure (measured at the wall orifice) as 
recorded for each location of the static probe as it was moved across the boundary layer. 

It is apparent that the wall static pressures are affected by the static-pressure probe when 
it is as far from the wall as 3,8 centimeters (1.5 in.), whereas the pitot-probe interference 
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had practically disappeared when the latter probe was 2.29 centimeters (0.9 in.) from 
the wall (fig. 2). As mentioned previously, these disturbed values of wall static pres- 
sures have not been used in any data-reduction procedures. 

Static pressures measured across the boundary layer with a static -pressure probe 
at x = 2.083 meters (82 in.) are also shown in figure 7(a). The nominal stagnation con- 
ditions for these data are shown in the figure. A sketch of the static probe and its location 
with respect to the nozzle-wall orifice is shown in the figure. Four orifices were located 
5.083 centimeters (2 in.) downstream of the tip and 90° apart around the probe circum- 
ference to reduce angle -of-attack effects. The accuracy and reliability of the data may 
be severely limited because of the large Mach number gradients in the boundary layer. 

The resulting change in Mach number across the 0.3 180 -centimeter -diameter (1/8-in.) 
probe was approximately 0.75. Nevertheless, the measured static-pressure data were 
corrected for bluntness and viscous induced pressure effects by the method of refer- 
ence 50 (eq. (49) from ref. 50) . No transverse curvature effects were included in the 
two-dimensional correction of reference 50. As a first approximation, the local Mach 
numbers used in this method were computed from the corrected pitot pressures (fig. 4) 
and the computed free -stream static pressure since the viscous induced corrections would 
be the largest in the outer part of the boundary layer where local Mach numbers are large. 
The local values of T t (from measurements corrected for radiation and conduction 
losses) through the boundary layer were also used. The surface temperature of the static - 
pressure probe was estimated by assuming a linear variation with y between assumed 
values for the ratio of probe temperature to free -stream total temperature of 0.2 at y = 0 
and 0.6 at y = 9.53 centimeters (3.75 in.). 

The resulting corrected static -pressure data shown in figure 7(a) indicate a trend of 
increasing pressure as the wall is approached. The same general trend is also indicated 
by the inviscid calculation but the latter distribution is displaced outwards from the wall 
by about 5.08 centimeters (2 in.). Also shown for comparison is a viscous correction to 
the measured data utilizing the first-order term of an empirical expression obtained from 
reference 51. The second-order term of this reference was obtained for a static -pressure 
probe similar to the present one but with a maximum y value of only 1.4. Hence, this 
second-order correction was not used herein since the present y values were about 14. 
The important point here is that the corrected static pressures at 6 (y ~ 9.53 centimeters 
(3.75 in.)) are in reasonable agreement with the value computed from the measured 
Pt 2 /P 0 ’ anc * a ^ so with the values from the inviscid -flow calculation. It is concluded on 
the basis of these results that there is a significant increase in static pressure from the 
edge of the boundary layer to the wall. 

The exact distribution of static pressure across the boundary layer cannot be 
determined from these data, because of possible wall -interference effects for 
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y % 3.81 centimeters (1.5 in.) and unknown effects of the Mach number gradient across 
the boundary layer. Another source of uncertainty is the thermal -diffusion effect on the 
static -pressure measurement. Estimates by the method of reference 36 indicate that the 
thermal -diffusion effect would increase the measured pressures by at most 3 percent for 
the range of probe wall temperatures utilized above in the viscous interaction correction. 
Because of these uncertainties, linear distributions between the measured wall value and 
the value inferred at y = 6 from measured p t 2 /P 0 have been used to obtain all of the 
following computed boundary -layer profiles. (The wall and stream static pressures used 
are given in table I.) This distribution is specified by the relations 
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Values of A have been limited to 1.0 and 0.28 which give, respectively, a simple linear 
distribution from y = 0 to 6, and a two-step linear distribution that approximates the 
data somewhat better. These distributions are compared with the corrected data in fig- 
ure 7(b) where p/p e is plotted against y/6. Shown for comparison are static -pressure 
measurements across the boundary layer on the wall of a helium tunnel at Mach num- 
ber 21.6 (ref. 11). The same trend of increasing static pressure as the wall is approached 
is evident in these data. 


Mach number distribu tions.- The local values of Mach number as computed from 
the Rayleigh pitot formula with ratios of local pitot pressure to total pressure and ratios 
of static pressure to total pressure are plotted against y in figure 8. Figure 8(a) shows 
the Mach number distribution out to about 18 centimeters (7 in.) from the wall as based 
on corrected pitot pressures (fig. 4) and static pressures from equation (1) with A = 1.0 
for all three conditions and with A = 0.28 for the intermediate condition 2 where 
p Q * 4310 N/cm^ (6250 lb/in^). The maximum effect of these two static -pressure dis- 
tributions was to change the Mach number by about 1.0 in the central region of the bound- 
ary layer. However, because of the previously mentioned uncertainties in the static - 
pressure distributions, the two-step linear distribution has been used only with the 
intermediate condition 2 for the purpose of illustrating the possible effects of different 
static -pressure distributions on the various profile parameters. The Mach number 
distributions tend to be nearly linear for y > 2.54 centimeters (1 in.) and the slope 
decreases with decreasing total pressure. 

Shown in figure 8(b) are the local values of Mach number in the wall region 
(y < 1.27 centimeters (0.5 in.)) for condition 2 (table I) as obtained with both the cor- 
rected and uncorrected pitot pressures. The solid faired line is used to evaluate 
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(a) Distribution across boundary layer and adjacent free stream. 

Figure 8.- Mach number distributions across the boundary layer. Equation (l) used for static pressures 






(b) Variation of local Mach number in the near-wall region from both uncorrected and 
corrected pitot pressures. A = 1.0. Flagged symbols show corrected values. 

Figure 8.- Continued. 

experimental skin friction as indicated in a subsequent section of this report. This line 
was faired through all data points and therefore represents the best estimate of Mach 
number distribution if equal "weight” or reliability is assigned to all data points. Com- 
parison of Mach numbers from corrected and uncorrected pitot pressures shows, how- 
ever, that the rarefaction effects cause larger percentage corrections in Mach number as 
the wall is approached. The skin friction (which is proportional to 9M/9y) would there- 
fore be strongly influenced by uncertainties in both static pressures and pitot pressures. 


34 



Predicted Mach number distributions from the finite -difference method of refer- 
ence 5 are also shown in figure 8(b). These theoretical results will be discussed in detail 
in a subsequent section and are included here to illustrate the significant differences 
between predictions of skin friction for laminar and turbulent boundary layers. 

In order to obtain some indication of the validity of the pitot -pressure corrections, 
some additional pitot-pressure surveys were made in the outer part of the boundary layer 
for y > 5.08 centimeters (2.0 in.) with large circular tubes of 0.3 18 -centimeter 
(0.125-in.) diameter. The results are shown in figure 8(c) where the local Mach num- 
bers computed from the pitot pressures measured with the large probe and the flattened 
probe are shown. Comparison of corrected and uncorrected data in the overlap region 
of 7.6 > y > 5. 1 centimeters (3.0 > y > 2.0 in.) indicates that while the corrections are 
generally of opposite sign for the two different tubes, the net scatter in the data is 
reduced. This reduction in data scatter indicates that the corrections are satisfactory 
for the limited range of conditions within the overlap region. Also, comparison of the 
corrected data from the large circular probe and the small flattened probe indicates that 
there was no significant effect of the Mach number gradient on the pitot pressure from 
the large probe for y > 5.6 centimeters (2.2 in.). Since the Mach number gradient is 
nearly linear, it can be concluded that the large probes could have been used over most 
of the boundary layer in the present investigation, however, all data presented in the 
tables and other figures are from the small flattened probe. 

Total -temperature surveys . - The ratios of local stagnation temperature to mea- 
sured settling -chamber temperatures are given in figure 9. Figures 9(a) and 9(b) show, 
respectively, the uncorrected and corrected data plotted against y. (See appendix A for 
details of the method used to obtain total -temperature corrections.) The nominal test 
conditions shown in figure 9 were taken from the temperature -survey tests (table I). 

These data are the first published experimental measurements, known to the present 
authors, of stagnation-temperature profiles in cold-wall turbulent boundary layers for 
Mach numbers as high as 19 in air or nitrogen. Comparisons between theoretical calcu- 
lations and some of the present data (without pitot -pres sure and final temperature cor- 
rections) have previously been made in references 5 and 47. The only other high Mach 
number (M % 19) data published to date have been in helium (refs. 11, 16, 26, and 27). 

The thermal boundary -layer thicknesses for the present data (shown in fig. 9(a)) 
are somewhat larger than the thicknesses based on pitot -pressure data (see fig. 4) also 
shown in the figure. This same result for relative thermal and pitot -pressure boundary - 
layer thicknesses has been observed in other experiments (private communication from 
the Von Karman Gas Dynamics Facility, Arnold Engineering Development Center). The 
wall temperature measured near x = 2.083 meters (82 in.) where the nozzle-wall sec- 
tion was water cooled was about 300 K (540° R) which gives a ratio of wall -to -stream 
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stagnation temperature of about 0. 17. No upstream wall temperatures were measured 
for the present tests. However, the wall temperatures probably do not exceed 330 K 
(600° R) for the water-cooled section of the nozzle (which extends from the end of the 
throat section at x = 14.5 centimeters (5.7 in.) (ref. 30) to the exit of the nozzle) since 
the heating rates are small downstream of the throat section. 

Comparisons of the uncorrected and corrected data indicate that the corrections 
caused large changes in both the magnitude and profile shapes of total temperature. 
Comparisons of theoretical predictions from the methods of references 5 and 52 with the 
corrected data (fig. 9(b)) again indicate that the boundary layer is not laminar. Further 
implications of the theoretical predictions will be discussed in a subsequent section of 
the report. 

The corrected total -temperature data are plotted as a function of yj 6-p in fig- 
ure 9(c) to show the consistency in profile shapes of these data. The "hump" in the pro- 
file data at y/6-p ~0.3 is present in all the corrected data and may be caused by or 
associated with the thick viscous sublayer which extends out to y/8-p - 0.2 as indicated 
in the figure. Since the unc -rrected data (fig. 9(a)) show little, if any, evidence of such 
a hump, it is apparent that the nonlinear effects of heat conduction and radiation losses 
on the indicated probe temperatures have smoothed out and removed the hump. 


Computed Boundary -Layer Profiles and Parameters 

The velocity profiles through the boundary layer were obtained by using the Mach 
number (fig. 8(a)) and total -temperature profiles (fig. 9(b)). With the assumptions of an 
ideal gas (p = pRT) and constant ratio of specific heats, the velocity profiles were then 
calculated from the equation 


_u_ = M i/T~ 
u e i/T e 
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and y = 7/5. The Mach number values used to calculate the velocities were taken from 
fairings of the data in figure 8(a) at values of y corresponding to the total -temperature 
data points of figure 9. The local mass flow was calculated from the relation 
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where p was computed from equation (1). 
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The displacement -thickness factor that appears in the integral equations for an 
axisymmetric boundary layer is 
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Similarly, the bovindary -layer momentum, thickness is 


( 5 ) 



The values of 6* and 6 for the present tests were obtained by numerical integration 
of equations (5) and (6). 

The values of skin friction for the present tests were obtained from the Mach num- 
ber gradient at the wall by using the formula 
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Values of (aM/9y) w were obtained from fairings of Mach number data near the wall 
similar to the solid line fairing shown in figure 8(b) for condition 2 (table I). However, 
the values of c f obtained from such fairings must be considered as only approximate 
estimates because of the large pitot-pressure corrections required in the wall region, 
uncertainties in the static pressures, and the relatively large spacing between data points 
near the wall. 

The Reynolds analogy factor has also been estimated from the total temperature and 
Mach number gradients at the wall by the formula 
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and 



0.9 


was assumed. While the results obtained from equations (7) and (8) are considered pre- 
liminary due to the uncertainty in the experimental values of the total temperature and 
Mach number slopes at the wall, they are included here as a matter of interest and for 
comparison with previous data and theoretical predictions. The numerically integrated 
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values of 5*/ 6 and d/6 from equations (5) and (6) and of Cj and 2Ngj-/cf obtained 
from equations (7) and (8) are given in table V(a). 


COMPARISONS OF EXPERIMENTAL RESULTS WITH 
THEORETICAL PREDICTIONS 

The corrected pitot -pressure and total -temperature data have been used to compute 
the profile parameters u/u e , M/M e , and 0 and the results may be compared with 
theoretical predictions in figure 10. Predicted values of M and T^/T 0 are also shown 
in figures 8(b) and 9(b) and predicted values of boundary -layer -thickness parameters, skin 
friction, and Reynolds analogy factors are given in table V(b) where the corresponding 
experimental values at the intermediate pressure conditions are given for comparison. 

As previously discussed, the local static pressures computed from equation (1) have been 
used to calculate the profile data points with A = 1.0 except for the intermediate stag- 
nation pressure run where A = 0.28 was also used. Before further discussion of com- 
parisons between experimental results and theoretical predictions, the principal ingre- 
dients of the theoretical calculations will be presented. 

TABLE V.- BOUNDARY-LAYER PARAMETERS 


(a) Experimental data 
[x = 2.083 meters (82 in.)] 


M e 

Average 

N Re/ m 

Average 

N Re/ft 

6 

676 

e/5 

6*/e 

(9M/ay) w 

c f 

2N St/ c f 

e 

N Re, 9 

cm 

in. 

per m 

per ft 

cm 

in. 

19.28 

1.96 x 10 6 

5.98 X 10 5 | 

11.2 

4.4 

0.530 

0.0159 

33.4 

295 

90 

9.0 x 10" 4 

0.97 

0.178 

0.070 

3490 

19.42 

2.76 

8.41 

9.4 

3.7 

.492 

.0185 

26.7 

240 

73 

5.2 

1.09 

.173 

.068 

4760 

19.65 

2.99 

9.12 

8.65 

3.4 

! 

.499 

.0179 

27.8 

243 

74 

4.3 

1.06 

.155 

.061 

4640 


(b) Comparison of theoretical predictions for turbulent flow 
with experimental values 

k = 19.48; N Re /m = 2.69 x 10^ (N Re /ft = 8.2 x 10 5 )] 


Source 

5 at 

M = 0.995 

N Re,6 

N Re,6 

e 

c f 

— 

cm 

in. 

cm 

in. 

Bushnell and Beckwith 

18.10 

7.14 

4.9 X 10 5 

7.8 X 10 3 

0.29 

0.114 

5.2 X 10" 4 

(ref. 5) 

Harris (ref. 52) 

9.53 

3.75 

2.57 

2.70 

.098 

.0385 

4.64 

Experiment 

9.40 

3.7 

i 

2.6 

4.8 

.173 

.068 

5.2 
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(b) Mach number. 


Figure 10.- Continued. 






M e Condition A Ref. 

(.Table 1) 



(c) Total -temperature-velocity profile. 
Figure 10.- Concluded. 


Theoretical Methods and Inputs 

The theoretical results shown in figures 9(b) and 10 were calculated for compress- 
ible laminar and turbulent boundary layers by two independent theoretical methods. The 
method of Bushnell and Beckwith (ref. 5) employs a two -point implicit finite -difference 
procedure. The method of Harris (ref. 52) uses a three-point implicit finite -difference 
procedure. Both methods solve the mean flow boundary -layer equations and utilize eddy- 
viscosity and eddy -conductivity (as determined from specified turbulent Prandtl number 
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distributions) concepts for the turbulent flow solutions. A simple mixing-length approach 
is used in the method of reference 5. The eddy -viscosity expressions used by Harris 
(ref. 52) are identical to those of Cebeci and Smith (ref. 53) who used the Prandtl mixing- 
length relation near the wall and a Clauser eddy -viscosity expression in the far -wall 
region. The calculation by the method of reference 5 incorporated the complete pressure- 
gradient history in the entire nozzle with appropriate initial profiles upstream of the 
throat, whereas the calculation by Harris was for an equivalent flat plate without pres- 
sure gradient. Comparison of the results from the equivalent -flat -plate solution with 
those from the solutions by the method of reference 5 will then provide some indication 
of the effects of the upstream history of the nozzle flow on local profile shapes near 
the nozzle exit. The nominal inputs for the flat -plate calculation were M e = 19.5, 
p e = 11.0 N/m^ (0.0016 lb/in^), and T e = 22 K (40° R). The inputs used for the method 
of reference 5 are given in figure 11. The wall temperature is the same as that of refer- 
ence 5, but the final Mach number at the 2.083-meter (82-in.) station has been increased 
from 18.2 to 19.5 in accordance with more recent experimental data. The stagnation 
conditions used in all solutions were p Q = 4310 N/cm2 (6250 psia) and T 0 = 1683 K 
(3030° R). 


Discussion of Comparisons 


The theoretical results calculated by the method of Harris (ref. 52) were selected 
from the equivalent flat-plate solution at the streamwise location where Np ej g (and 
hence 6) agreed with experiment. The computed value of Np e> q would then agree with 
the experimental value only if both the Mach number and velocity profiles agreed perfectly 
with data. Figures 10(a) and 10(b) show that such agreement was not obtained; hence, 
agreement with the experimental value of Np^g was not obtained. (See table V(b).) 

The good agreement between predicted skin friction from this solution and the experi- 
mental value given in table V(b) is probably fortuitous since the values of Np e g are not 
in agreement. 


The results from the method of reference 5 (figs. 9(b) and 10) were taken from the 
nozzle calculation at x - 2.083 meters (82 in.) corresponding to the experimental sur- 
vey station. The predicted boundary -layer thickness from the calculation for turbulent 
flow is almost twice as large as the experimental value, whereas the predicted value of 
6 for laminar flow is about 2/3 of the experimental value. The calculated values for 
turbulent flow of Q, Np e ^g, and Np e , g are therefore too large (see table V(b)) and 
again the agreement in skin friction is fortuitous. ^The theoretical skin-friction values 
for all solutions were computed with the approximation that ftlV/r M 


/ 9M\ ,, M 

v 9 y/ w y 


evaluated from 


was 


the solution at y ~ 0.13 centimeter (0.05 in.). This evaluation for (aM/9y) w 
intended to approximate the experimental evaluation which is based essentially on the 
two or three data points nearest the wall as illustrated in fig. 8(b).j The solutions for 
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this same nozzle flow as reported in reference 5 gave boundary -layer thicknesses from 
9.9 to 12.7 centimeters (3.9 to 5.0 in.) depending on the modifications to the eddy vis- 
cosity and turbulent Prandtl number functions as well as the input and boundary condi- 
tions used. The eddy viscosity and Np r T functions used to obtain the profiles shown 
in figures 9(b) and 10 are identical to those for solution number 5 of reference 5 for which 
the predicted value of 6 was 12.7 centimeters (5.0 in.). However, the boundary con- 
ditions for the present solution were modified downstream of the x = 1.22-meter (48-in.) 
station, to provide inviscid-flow conditions that are consistent with stagnation conditions 2 
(table I) and with the present measurements of stream static pressure. The resulting 
inputs are compared with values used in the corresponding solutions of reference 5 in the 
following table: 



P, 


Pol/ 2 c p^o 

x = 2.083 meters (82 in.) 

Source 

D 

Mo 

p'p n 

M e 

t 

) 


N/cm 2 

lb/in 2 

per cm 

per in. 

B( O 


cm 

in. 

Previous data 
(solution 5, 
ref. 5) 

4850 

7040 

3.15 x 10 7 

0.80 x 10 8 

2.47 x 10“ 5 

18.2 

12.7 

4.98 

Present data 
(figs. 9(b) 
and 10) 

4310 

i 

6250 

2.80 

.71 

! 

1.93 

19.5 

18.1 

7.14 


Subsequent comparisons of the two solutions have shown that downstream of the 
x = 1.2-meter (48-in.) station the boundary -layer thickness varied inversely as the 
density ratio while the length of flow had little effect on the rate of boundary -layer 
growth, at least for flow distances of several boundary -layer thicknesses. These results 
indicated that the local density must be accurately known before reliable predictions of 
boundary -layer thicknesses can be obtained for hypersonic nozzles. 

In figure 10(c), the normalized total -temperature parameter 0 is plotted as a 
function of the velocity ratio u/Ug. The prediction by the method of reference 5 is in 
good agreement with the data for the outer part of the boundary layer (u/u e > 0.8). The 
flat-plate solution is nearly the same as the linear variation of Crocco which is in good 
agreement with the data for u/u e < 0.8. The differences in these predictions are caused 
partly by the upstream pressure gradient history which is incorporated in the nozzle flow 
calculation but not in the flat -plate calculation. The magnitude and distribution of the 
turbulent Prandtl number also affects the variation of ~0 with u/u e , particularly for 
u/u e < 0.8, as shown in reference 5. 
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Comparisons of predictions from both methods for turbulent and laminar boundary 
layers with the experimental velocity and Mach number profiles (figs. 10(a) and 10(b), 
respectively) indicates the boundary layer on the nozzle wall is not laminar at this 
x = 2.083-meter (82-in.) station. From the previous discussion of figure 10 and the 
results shown in reference 5 for this nozzle flow, it is evident that not only must the 
local density levels be accurately known, but also the entire flow history must be included 
in the calculation before agreement with experimental profiles and boundary -layer thick- 
ness can be expected. Nevertheless, the equivalent -flat -plate solution for turbulent flow 
gives much better agreement with Tt/T 0 , u/u e , and m/m 6 profiles (figs. 9(b), 10(a), 
and 10(b), respectively) than the solution with the upstream flow history included. Of 
course, the boundary -layer thickness is forced to agree by matching Npj e § with the 
data so that no prediction of 6 is obtained. Other than this lack of any prediction for 
boundary -layer thickness, the only notable deficiency in the flat -plate solution is the 
overprediction in T^/Tq for 3.8 cm < y < 7.6 cm (1.5 in. < y < 3 in.) (fig. 9(b)) and 
the resulting overprediction of ~9 for u/u e > 0.8 as shown in figure 10(c). The nearly 
quadratic trend of 9 with u/u e indicated by the data in this region are characteristic 
of nozzle -wall boundary layers (ref. 47) and are believed to be caused mainly by upstream 
pressure and wall temperature gradients. The nearly linear trend of the present data 
for u/u e < 0.8 is typical of flat -plate data and may be related to the thick viscous sub- 
layer which, from figure 10(a), extends out to values of u/u e = 0.7 to 0.8. 

Since the predictions for the profile shapes by the turbulent flat -plate solution 
are in reasonable agreement with data, the most likely cause of the overprediction in 
boundary -layer thickness at x = 2.083 meters (82 in.) from the solution by the method 
of reference 5 (obtained with the proper upstream boundary conditions) would be simply 
an overprediction of the rate of boundary -layer growth in the upstream region of the noz- 
zle. One way to obtain predictions of smaller growth rates is to reduce the eddy viscos- 
ity by large amounts. The effects of various modifications to the eddy viscosity and other 
inputs used in the theory will now be considered in an attempt to identify the cause of the 
overpredictions in boundary -layer thickness. 

Effects of Modifications to Theory Inputs 

The eddy viscosity is defined in terms of the turbulent shear and normal velocity 
gradient as 


t T = “(P v )' u ' = e 


9u 

ay 


( 9 ) 
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Expansion of the Reynolds stress correlation 


(pv)'u’ = p v'u' +vp'u' + p’v’u* 

and introduction of the mixing length l and intermittency functions in the x- and 
y-directions, r x (x) and r y (y), results in the following expression for e (see ref. 5): 


e - r x r y 


p; 2 

au 


9U 

)1 


ay 

9y\ 

9y 

/J 


( 10 ) 


where 


6 


1 - exp 


yp w 


'w 


: ( x > I) 


A*Pw rw. 

The mixing-length function for the near -wall region is 
fn = 0.4 \ 

and for the far -wall region 

ff = Constant = 0.09 

An alternate form used in reference 5 for ff is 
f f = f f (Hi*) = 0.1234Hi* - 0.0317 


(ID 


}S0.1 

0 


?■ = 0.3 

o 


( 12 ) 


(13) 


\ ^ 0.3 (14) 

o 


which was designated "variation 2” in reference 5. The values of f for the range in 
y/6 of 0.1 £ y/6 S 0.3 were determined from a straight-line segment between the 
limits for f n and ff as noted previously. 

Laminar solutions are obtained for e = 0 or T x (x) = 0. If an established turbu- 
lent boundary layer is subjected to large streamwise accelerations, the turbulence may 
be suppressed, and T x would be correspondingly reduced. Several criteria for such 
laminarization effects have been suggested, one of which is 


K = 


_ Me du e 
P e u e 2 dx 


In the investigation of reference 54, when K exceeded about 2 to 3 x 10 _6 in the approach 
to the nozzle throat, the boundary -layer profiles and heat-transfer rates in the throat 
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region tended to approach those for a laminar boundary layer. This laminarization 
trend persisted downstream of the throat minimum, even though the local values of K 
decreased below the level of 2 x 10 “ 6 . 

The variations with axial distance of K, Nr 6jX , and Njj e> q (for laminar and 
turbulent solutions) for the present nozzle are given in figure 12. It is seen that K 
exceeds the "laminarization criteria" upstream of the throat and hence the wall bound- 
ary layer presumably could be partially or completely laminar in this region. If the 
boundary layer did revert to laminar flow, this condition might persist far downstream 
because of the associated small values of Nr 6? q which peak at 3000 in the throat and 
decrease to a minimum of 1400 about 30.5 centimeters (12 in.) downstream of the 
throat. On the basis of these levels and distributions of K and Nr 6 q a trial T x 
distribution shown in figure 11 was used for several solutions. 

The Ty(y) function in equation (10) is used to represent the effect of intermittency 
in the outer part of a turbulent boundary layer. The density gradient term in equation (10), 
referred to herein as the p' term, decreases the eddy viscosity in the outer part of the 
boundary layer if 9p/9y is positive. The eddy viscosity is also decreased by the use 
of either of the two intermittency functions, which vary between 0 and 1.0, and by reduc- 
tions in mixing length. 

Several modifications to the eddy viscosity of this nature were used in various com- 
binations to obtain a series of solutions for the present nozzle flow. In order to illustrate 
the effect of such modifications, results from four typical solutions are presented and 
compared with experimental results in figure 13 and table VI. Solution number 1 is the 
sa’'’e as the turbulent solution already presented in figures 9(b) and 10 and in table V(b) 
by he method of reference 5. The turbulent Prandtl number distributions used in these 
four solutions are given in the inset in figure 13(a) and other inputs are listed in table VI. 
Nore of the four solutions shown utilized the p' term because it was found that while 
the eddy viscosity and boundary -layer thickness were reduced the profile shapes were 
not sufficiently full. 

TABLE VI.- EFFECTS OF EDDY VISCOSITY MODIFICATIONS AND COMPARISON WITH EXPERIMENTAL RESULTS 








Values at x = 2.083 meters (82 in.) 




f f 

r x 

r y 

(Z/6) m ax 

-Hi* 

6 


c f 


9 

6 

* 

6*/0 

N Re,0 




cm 

in. 

cm 

in. 

cm 

in. 

Equation (14) 

1.0 

1.0 

0.069 

1.724 

18.1 

7.14 

5.22 X 10" 4 

0.29 

0.114 

12.20 

4.807 

42.1 

7810 

Equation (14) 

Figure 11 

1.0 

.055 

2. 143 

9.87 

3.89 

4.43 

.130 

.051 

7.15 

2.818 

55.5 

3470 

0.09 

Figure 11 

1.0 

.09 

1.918 

11.65 

4.59 

5.56 

.178 

.070 

7.18 

2.823 

41.6 

4770 

.09 

Figure 11 

Reference 5 

.09 

1.920 

10.70 

4.21 

5.29 

.160 

.063 

7.86 

3.098 

49.2 

4310 






9.4 

3.7 

5.2 

.173 

.068 

4.63 

1.82 

26.7 

4760 






8.64 

3.4 

4.3 

.155 

.061 

4.31 

1.70 

27.8 

4640 
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(d.) Total-temperature parameter variation with velocity. 
Figure 13.- Concluded. 


The only way to reduce the boundary -layer thicknesses significantly while retaining 
good profile shapes was by using r x = 0 over some upstream portion of the nozzle. 

The Np r p function has little effect on boundary -layer thickness but does control the 
Crocco profile directly. These effects of the T x function and Np r p are illustrated 
by comparisons of results from solutions 1 and 2. While solution 2 has about the right 
boundary -layer thickness, the total -temperature and velocity profiles (figs. 13(a) 
and 13(b)) are not sufficiently full for y < 5.10 centimeters (2 in.). The Mach number 
and Crocco profiles (figs. 13(c) and 13(d)) are considerably improved. In order to obtain 
fuller profiles for y < 5.10 centimeters (2 in.), (Z/6) max was increased by using 
ff = 0.09, and further reductions in Np r T gave good agreement with the Crocco profile. 
Comparison of theoretical profiles and 6, Cp and 6 from solutions 2, 3, and 4 with the 
experimental results (fig. 13 and table VI) shows that these solutions essentially bracket 
the data except for the Mach number profiles. The large range in eddy viscosity for these 
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Figure l4.- Eddy-viscosity distribution at x = 2.083 meters (82 in.) 

theoretical solutions of figure 13. 


from the four 


four solutions is shown in figure 14. These distributions should not be used to infer any 
general causative trends between eddy viscosity and flow variables since the eddy- viscosity 
values were obtained from the solutions at x = 2.083 meters (82 in.) and are influenced 
by the entire upstream history of the boundary layer as well as by local conditions. 

The results shown in figure 13 and the preceding discussion indicate that further 
adjustment of the theoretical inputs, particularly T x , Np rj -p, and (z/6) max , the pre- 
dictions could probably be brought into as close agreement with the data as desired. How- 
ever, further adjustments and solutions are not justified because of limitations in the data 
and theory. 
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The principal limitations in the data, besides the usual experimental uncertainties, 
are the unknown extent of upstream laminar ization and the unknown wall -temperature 
history. Comparison of the preliminary version of the present data as given in refer- 
ence 5 with the final corrected data herein emphasizes the experimental uncertainties 
that are of particular concern for all boundary -layer data obtained at hypersonic condi- 
tions. As mentioned previously, the present data as published in reference 5 did not 
include the rarefaction corrections to pitot pressure, and the total -temperature correc- 
tions were preliminary estimates. Also, the previous free -stream density was almost 
30 percent larger than the present value for a stagnation pressure roughly 10 percent 
smaller. This reduction in free -stream density for the present conditions increased the 
predicted boundary -layer thickness by about 40 percent. Thus, all of the present con- 
clusions concerning the comparisons between predictions and experimental results hinge 
on the accuracy of these large corrections. It is believed that the "final" corrected data, 
as presented herein, are accurate to within ±10 percent. The new corrections have 
changed the preliminary results as presented in reference 5 only for y/5 < 0.5 where 
the maximum changes were: an increase in Tt/T 0 of about 20 percent, a decrease in 
M/M e of nearly 50 percent, and a decrease in velocity of about 20 percent. These 
changes in Tt/T 0 and u/u e produce a nearly linear variation in 9 with u/u e (for 
u/u e < 0.8) rather than the quadratic type variation shown previously in reference 5. The 
quadratic variation typical of most previous nozzle-wall boundary layers may be caused 
partly by cold-wall conditions in the settling-chamber and nozzle-throat regions of most 
facilities, a condition which is not present in the Mach 19 nitrogen facility. Recent data 
obtained by W. V. Feller in the Langley Mach 6 facility show that the typical quadratic 
Crocco profile near the nozzle exit is modified to a more nearly linear profile by 
increases in wall temperature in the settling chamber and throat region of the nozzle. 

The principal limitations in the theory are the same as those for any turbulent 
boundary theory, namely, the unknown relations between turbulent flux correlations and 
mean flow quantities. At present, the best that can be done is to model these correlations 
and essentially by trial-and-error processes, as demonstrated herein, to determine what 
modifications to the models give the best predictions. The real problem with this process 
is that the results may not be unique. Another limitation in the theory is the assumption 
of constant static pressure which is believed to be the main reason for the overpredictions 
in 6* (see table VI) from solutions 2 to 4. The overpredictions in 6* may also be 
partly caused by transverse curvature effects, which were included in the theory, but only 
to first order. (See ref. 5.) Nevertheless, on the basis of these comparisons, it seems 
likely that the boundary layer was laminarized to some extent in the throat region. It is 
also apparent from the present results and those of reference 5 that in order to obtain 
agreement with the present experimental Crocco profiles, values of Np r p of 0.5 or 
less are required over the inner part of the boundary layer. Again, due to uncertainties 
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in the data and in the upstream boundary conditions and r x variation, these conclusions 
must be considered as only tentative. 

COMPARISONS WITH PREVIOUS DATA 


To the authors 1 knowledge, corrections for viscous or rarefaction effects on pitot- 
pressure measurements have not been applied to the previous data used in the following 
comparisons except the data of reference 11. 

Mach Number Profiles 

The present Mach number profiles are compared with previous data for nozzle -wall 
boundary layers in figure 15. The data in this figure illustrate the large differences that 
can occur in profile shapes over the wide range of test conditions represented. It can be 
seen that, in general, the profiles become less full as the Mach number is increased and 
N Re , q is decreased. In fact, it was found that the parameter M e j ^Np^T g generally 
increased as the profile shapes shown decreased in fullness. Thus, consideration of the 
data for y = 7/5 shows that the fullest profiles (those to the left in the figure, refs. 12 
and 28) have 0.03 ~ M e /|/Njj e g x 10^ ~ 0.1, while the group of profiles that are more 
nearly linear (refs. 15, 17, and the present data) have 0.2 ~ M e /|/Np^ e g x 10^ ~ 0.3. 

The data of Wallace (ref. 19) have M e /^Npj e g values of about 0.11 to 0.14 which are 
between the aforementioned ranges, and the profile shapes are consistent with this range 
since they follow the profile group with full shapes for y/5 <0.1 and change over to the 
"linear" group for y/5 > 0.4, The helium profiles follow essentially the same trends 
as noted for the air and nitrogen data. 


Velocity Profiles 

The values of N in the power -law velocity profile relation 



(15) 


are of use in characterizing and classifying velocity profiles. 

The values of N which best match the outer portion of the present velocity profiles 
are shown in figure 16. These values of N have been compared with results from pre- 
vious data for axisymmetric nozzles in reference 55. The present data for N are near 
the upper side of the data band shown in reference 55 where it was indicated that for 
Nr g , q ~ 2500, N would increase approximately as ( T w/ To) 0 * 25 . In order to determine 
if this variation is also present at other values of Nr 6 q, values of N(T w /T q ) D from 
the present and previous data (data for adiabatic wall temperatures and M e > 17 have 
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Figure 15.- Comparison of experimental Mach number profiles with previous results on nozzle 
walls for 8.25 ^ Me = 37 and 1.9 X 10? < NR e ,9 < X 1CP. Data from reference 11 
are the same as data of reference 29, hut corrected for pitot-pressure errors, and vari- 
able local static pressure. 
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TABLE VH.- REFERENCES AND VALUES OF EXPERIMENTAL PARAMETERS USED IN FIGURES 17 AND 19 


Reference 

y 

M e 

N Re,6 

N Re, 9 

T w/ T e 

T w 

/ T o 

(®s/ 6 )meas 

N 

O 15 

7/5 

11.85 

3.12 x 10 4 

0.94 x 10 3 

9.10 

i 0.312 

0.350 

6.45 


11.99 

4.21 

1.26 

9.89 


332 

.325 

6.50 



12.07 

5.25 

1.63 

10.40 


344 

.300 

6.57 

□ 14 

7/5 

6.65 

2.25 

2.79 

3.12 


312 

.100 

3.31 



11.30 

4.49 

3.77 

7.27 


317 

.200 

3.69 



11.50 

4.70 

2.58 

7.36 


268 

.220 

4.92 

O 17 

7/5 

8.25 

18.20 

12.0 

1.93 


150 

.010 

4.10 

16.30 

7.88 

1.85 

5.58 


118 

.080 

4.16 



14.00 

6.98 

2.40 

4.07 


119 

.075 

4.09 



12.65 

6.55 

3.17 

3.25 


117 

.070 

2.71 

A 12. 

7/5 

8.87 

40.60 

16.7 

5.90 

.352 

.030 

6.00 


9.00 

79.00 

31.3 

5.43 

.312 

.025 

6.10 



11.27 

25.00 

6.63 

8.42 

.319 

.060 

6.70 



11.60 

40.40 

13.84 
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*Data obtained from the Von Karman Gas Dynamics Facility, Arnold Engineering Development Center. 



been excluded, see table VII) were plotted against Npj e g in figures 17(a) and 17(b). 
Comparison of these two figures show that the correlation of N over the complete 
range of g is somewhat improved (particularly for the present data) with D ~ 0. 1 

(fig. 17(b)) rather than D ~ -0.25 (fig. 17(a)) as would be indicated by the results of 
reference 55 at Nj^ e g ~ 2500. 
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(a) Effect of T w /T 0 on correlation of exponent N in power-law velocity profile 

according to reference 55. 


Figure 17- - Correlation of exponent in power-law velocity profiles. (See table VII 
for symbol key.) y = 7/5, 0.07 $ T w /T 0 = 1.0. 
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In figure 18 the present velocity -profile data are compared with data for the com- 
plete profiles from previous investigations of nozzle -wall boundary layers. This figure 
illustrates how the viscous sublayer thickness 5s (relative to 6) tends to increase with 
increasing Mach number, decreasing Njj e q and decreasing T W /T 0 . In order to deter- 
mine quantitatively this increase in 6 S , it was defined as the value of y at the edge of 
the region where the velocity profile is approximately linear, In figure 19, the resulting 
values of 5 s /5 are plotted against M e |j/Njj e> q , the same parameter that characterized 



(b) Improved correlation of present and previous data for velocity profile exponent N. 

Figure 17*- Concluded. 
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the shape of experimental Mach number profiles (fig. 15). This plot shows a trend for 
6 s /6 to decrease as M e j g decreases. There is apparently no consistent effect 
of T w / T 0 on the scatter of the data points about the mean variation shown. The values 
of these parameters for all data used in the figure are given in table VII. While the cor- 
relation of figure 19 shows considerable scatter, it should be useful for estimating the 
sublayer thickness which is important in relation to roughness size responsible for 
increases in friction drag. 
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Figure 18.- Comparison of present data with previous results. 
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Figure 1 9.- Variation of viscous sublayer thickness with M e ^Np, e q for nozzle-wall bound- 
ary layers. Symbol key, experimental parameters, and references given in table VII. 


Variation of Total Temperature With Velocity 

— Tt - T 

In figure 20 the total-temperature parameter 6 = — is plotted as a func- 

T t,e “ T w 

tion of u/u e from the present data and from previous data for a wide range of M e and 
N Re 9 • is seen that most previous data are considerably below the linear Crocco var- 
iation which best correlates flat-plate flows. (See ref. 47.) The nozzle-wall data shown 
in the figure are from some of the same sources used in reference 47 except the recent 
data from the Von Karman Gas Dynamics Facility, Arnold Engineering Development 
Center, and the Langley data of W. V. Feller which were not included in the survey of 
reference 47. The empirical relation H = (u/u e )2 seems to best correlate all previous 
nozzle -wall data except the new unpublished data obtained by Feller at Langley in the 
same Mach 6 facility used to obtain the data of reference 28. These new data were 
obtained with the settling-chamber walls and throat region heated to temperatures above 
the free -stream stagnation temperature. The results suggest that the quadratic variation 
of "0 with u/u e , considered typical of nozzle-wall boundary layers, is caused by the 
persistent effects on the developing boundary layer of the comparatively cold-wall tem- 
peratures in the settling chamber and throat. The arrangement of the heater element 
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within the settling chamber and the throat construction used in the present facility (see 
ref. 30) indicate that the effective surface temperatures in the settling chamber are some- 
what higher than stream stagnation temperature while the throat temperatures are prob- 
ably as high as 0.8 of total stream temperatures. Hence, the nearly linear variation of 
0 with u/u e for the present data are consistent with the new results of Feller in the 
limited sense that higher surface temperatures in the settling-chamber and throat regions 
cause a trend toward the linear variation. 


ESTIMATED VELOCITY FLUCTUATIONS 


Estimated values of longitudinal -velocity fluctuations have been made in refer- 
ence 45 for the data of reference 19 (M e ~ 8.5) and compared with previous measured 
values for 0 g M e = 4.7. The location of the peak in velocity fluctuations is near the 
viscous sublayer edge which tended to move away from the wall with increasing Mach 
number (ref. 45). An attempt has been made to estimate the variation in velocity fluctu- 
ations through the present boundary layer and the results are compared with those of 
reference 45. 


The first-order y-momentum equation for axisymmetric turbulent -boundary -layer 
flow may be expressed as 



cos 0 

r w - y cos /3 




(16) 


where 0 is the nozzle contour angle at x, and r c is the radius of curvature of the 
surface in the x -direction. At a given x -location, equation (16) becomes a linear, ordi- 
nary diffe rential equation for the variable pv’2. Hence, the solution of equation (16) 
for j/v'2 is 



This equation shows that in order to get real values of \j v’2, 3p/9y would have to be 

negative over most of the boundary layer since the second term in the integrand is always 
positive and generally small. 

If it is assumed that 6/r c and S/r w are very small corresponding to small cur- 
vature at the wall, equation (17) reduces to 
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( 18 ) 


which shows again that 9p/ay would have to be negative near the wall while elsewhere 
P < P w * 

Equations (17) and (18) were used to compute l/^2/ue and the results are shown 
in figure 21. Equation (1) was used for the static -pressure distribution in equations (17) 
and (18) with A = 1.0 and in equation (18) with A = 0.28. The values of /3, 6/r c , and 

6/r w are given in the figure and were based on the inviscid nozzle -wall coordinates at 
x = 2.083 meters (82 in.). The Mach number and density profiles used were those for 
condition 2 with A = 1.0. 


The present results with A = 1.0 are roughly similar in trends and levels to the 
incompressible data of Klebanoff (ref. 56) also shown in figure 21. The curvature terms 
have very little effect except near the outer edge of the boundary layer. This result indi- 
cates that in terms of the boundary -layer coordinate system used in equation (16), the 
static -pressure variation across the boundary is caused mainly by the turbulent velocity 
fluctuations normal to the surface. However, in the present tests, finite values of 8p/8y 
may be imposed upon the boundary layer by the external inviscid-flow field. (See fig. 7.) 
Hence, the relative contributions of turbulent fluctuations and the inviscid-flow field to 
the observed static-pressure variation cannot be determined. Any variation in static 
pressure caused by turbulent fluctuations is enhanced by large stream Mach numbers and 
large density changes across the layer as indicated by equation (18), and the data com- 
parisons of figure 7(b). 


The data of reference 56 indicate that the approximate relation 



(19) 


is applicable over m ost of the boundary layer except near the wall in the viscous sublayer 
region. Values of ]fu^/u e have been obtained by applying the \l2 factor t o the results in 
figure 21 for finite -curvature values. Another approximate relation for i/^2 

is based 

on Prandtl's mixing-length concept: 


(u^ _ l K u / u e) 
u e 6 9(y/5) 


( 20 ) 


This relation has also been applied to the present data for p Q = 4309 N/cm^ (6250 lb/in.^) 
with equations (12) and (13) used for 1/ 6. The velocity -profile derivatives were obtained 
from the faired line shown in figure 10(a). 
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Figure 21.- Comparison of estimated normal-velocity fluctuations with data of reference 56 . 

Present results calculated with p =3.5°. 


The resulting estimated values of u e from equations (19) and (20) are plotted 

as a function of y/6 in figure 22. Also shown for comparison in the figure are previous 
experimental values as presented in reference 45. Again, the present values from equa- 
tion (19) are roughly similar to previous results while the mixing -length estimates are 
larger except for y/6 > 0.5. Because of uncertainties in the present mean flow data 
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Figure 22.- Comparison of present estimates of velocity fluctuations with previous data. 

(particularly the static -pressure data), these results for fluctuating velocities are only 
indicative of possible trends and approximate order of magnitude levels for the present 
conditions. More reliable estimates for fluctuating quantities will have to be obtained 
from direct measurements of p', p, and improved measurements of static pressures 
and other mean flow quantities. 

CONCLUSIONS 

Total -temperature, pitot-pressure, and static -pres sure measurements were 
obtained in the boundary layer of a contoured nozzle at Mach 19.5 in nitrogen. The 
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ratio of wall temperature to settling-chamber temperature was about 0.17; the boundary 
layer was about 10 centimeters (4 in.) thick, and the Reynolds number based on momen- 
tum thickness varied from about 3500 to 4800. A general assessment and review of probe 
errors encountered in pressure and total -temperature measurements in hypersonic 
boundary layers is included. Accordingly, the pressure data were corrected for viscous 
and rarefied flow effects and the temperature data, obtained with a bare -wire thermo- 
couple, were corrected for heat losses by radiation and conduction. The significant con- 
clusions from this investigation were: 

1. A static -pressure gradient exists normal to the nozzle wall as indicated by mea- 
sured static pressures at the wall and in the boundary layer and calculated free -stream 
pressures. This static -pressure gradient was apparently caused partly by inviscid-flow 
expansions and partly by fluctuations in the normal -velocity component. 

2. Viscous and rarefied flow effects required large corrections to the pitot -pressure 
measurements as far as 5.1 centimeters (2 in.) from the wall. In general, the trend was 
to reduce the pitot pressure in the wall region since the rarefied flow effects increased 
the indicated pitot pressure above continuum values. 

3. Estimated values of laminar sublayer thickness and the exponent N in the 
power -law velocity profile were in good agreement with previous results. The results 
indicated the sublayer thickness generally increased with increasing Mach number and 
decreasing momentum -thickness Reynolds number N R e> 0 . The exponent N increases 
with increasing N Re> q and is affected only slightly by the ratio of wall temperature to 
settling -chamber temperature. 

4. The relation between the normalized total temperature and velocity was nearly 
linear over the inner part of the boundary layer as compared with a quadratic relation 
found in previous nozzle -wall boundary -layer data. This result was attributed to differ- 
ences in surface or near surface temperatures in the settling -chamber and throat regions 
of the present nozzle and typical nozzles of previous investigations. 

5. Comparisons of the present results in terms of boundary -layer -profile shapes, 
thickness parameters, and estimated values of skin friction with predictions from finite - 
difference theories showed that the boundary layer was not laminar. The theoretical 
results could be brought into reasonable agreement with the data only by suppressing the 
turbulent flow in the throat region of the nozzle and by using turbulent Prandtl numbers 
considerably less than unity in the theory. Although these theoretical results may not be 
unique because of nonlinear interactions between the various options used as inputs, the 
comparisons indicate that the nozzle boundary layer may have been laminarized in the 
throat region. 
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6. Estimates of the intensity of velocity fluctuations in the boundary layer were 
based on measured values of mean static pressure and also on an assumed mixing-length 
distribution. These estimates were roughly similar in trends and levels to previous 
experimental values for incompressible and supersonic flow. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 16, 1971. 
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APPENDIX A 


EFFECTS OF HEAT LOSSES BY CONDUCTION AND RADIATION 
ON BARE -WIRE THERMOCOUPLE PROBES 

By Ivan E. Beckwith, William D. Harvey, and Christine M. Darden 

Basic Equations 

The notation used in temperature data-reduction procedures for all test runs as 
reported herein is defined by the following sketch: 



x = -L/2 

x = L/2 



< 

X 

i> 

K- 

C Wire A 

() (alumel) 



Wire C Q L 

(chromel) ) , r 

| 

Ts,a 1 

1 t 

r M T S,C 


Sketch A 


where the temperatures Tg Tg q, and T^ (the junction temperature) have, in gen- 
eral, been measured. The heat -balance equation for an element of length dx of the wire 
is 


;s 5 * + - T)dx 4 P W C f dx - I f )dx + f (t* - |/L T a ^)dx (Al) 


The simplest way to evaluate the conduction term in this equation is to assume that the 
temperature distribution along the wire is a quadratic of the form 


T = A + Bx + Cx^ 


(A2) 


The conditions of continuity in temperature and heat conduction at the junction (x = 0) are 
then utilized to give the temperature distributions for the alumel and chromel wires, 
respectively, as 


T A = t M + b a x + ^( t s,a - t M + j b a ) x2 

(A3) 

T C - T m + B A Kx + ±(t S:C - T m - \ KB A ) X 2 

(A4) 
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APPENDIX A — Continued 


where 


and 


Application to Probe Corrections 

The local total temperature T^= ? T aw ) is then calculated from equations (Al), 
(A3), and (A4) by the following procedure: 

The conduction term in equation (Al) is first written as 



Evaluation of equation (Al) at the junction for wires A and C in sketch A then gives the 
relations 



IgA ^u,^ / n u £HM_!^A Ba 2 

nd 2 + d 2 l rl t 1 M j 4 (, p w c j A dt 4 dT b A 


^(r b a - TmFa ) - ( T « 4 - 


1 


(A5) 


and 


lec %u,t k t / ? _ \ _ i , , dT^ j dk c 2 2 

nd 2 d 2 l rT t Tm / 4l p w c ) c dt "4 dT KBA 


2k 


C,M L L 


+ TH K T : B A + T M F cj + -dH T M - T 


oe f 


where 


F A = 1 " ^ 

t m 


(A 6 ) 


F C = 1 - 

t m 


75 


APPENDIX A — Continued 


Subtraction of equation (A5) from equation (A6) then gives a quadratic equation for B A 
which is 


^( e C - e A) = f(Pc c C -PA^a) 


dT 


M 


dt 


+ 2 


5 A,M 


l/^A „2 ^ b 2 


4 \dT 


dT 


+ 


2T M 


k C,M F C 




(A7) 


This equation is solved for B A as a function of the measured junction temperature Tj^ 
and the measured end temperatures Tg A and Tg c which are shown in figure 23 in 
the form of the variables Tq/Tm and Tm/T 0 . These data were obtained with the 
probe shown schematically in the figure. This probe is the same size and shape as the 
boundary -layer probe shown in figure 1 but with 0.076-millimeter-diameter (0.003-in.) 
wires attached, as indicated, to obtain the end temperatures. The faired curves in fig- 
ure 23 for F a and F^ have teen used herein for the reduction of boundary -layer sur- 
vey data. All of these present data were obtained without electric heating and for time- 
steady conditions, so the terms containing I and dT-^dt in equations (A5) to (A7) 
were zero. The values for thermal conductivity of alumel and chromel were taken from 
the faired lines shown in figure 24 (from data in ref. 57). Shown for comparison in this 
figure are the values of thermal conductivity of similar "type K" thermocouple wires 
from reference 58. The emissivity difference in equation (A7) is taken from 

faired curves of measured emissivity as discussed in appendix B where it is shown that 
the emissivity of chromel tends to increase more than e A after similar exposure 
to high-temperature nitrogen. 

The values of B A obtained from equation (A7) are substituted into equation (A5) 
which is solved for Tj. by using data for r and NN Uj t as given in reference 59. 

These data apply only up to about Mach 5; however, preliminary measurements of the 
transient heating rates of typical boundary -layer probes during the present investigation 
indicate that the N^ u ^ curve of reference 59 for Mach 5 is still applicable at Mach 19. 
Since the emissivity of both alumel and chromel does depend on the surface condition of 
the wires, it is first necessary to determine a reference emissivity for each test. This 
value is obtained by solving equation (A5) for e A with all measured quantities at y = 5 
used and with T^- = Tj. e required. The resulting value of emissivity, termed e A e , 
is used to construct an effective emissivity curve for the remaining data from that partic- 
ular run. (For details of this procedure and the effective emissivity curves used, see 
appendix B.) 

To obtain T^/T^ e from equation (A5), an iterative procedure was required, since 
Nnu t an d f are functions of local Reynolds number, Mach number, and Knudsen num- 
ber (ref. 59). Since large corrections to the local Mach number were required due to the 
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Figure 23.- Variation of ratio of end temperature to junction temperature with junction temperature. 
Linear dimensions in inches (parenthetically in millimeters). 


APPENDIX A - Continued 


APPENDIX A - Concluded 


rarefied flow effects (see fig. 3 and discussion in text), a double iteration cycle was used 
to obtain the final corrected pitot pressures and total temperature. The accuracy criteria 
used in the iteration procedure was that the percentage change with successive iterations 
was less than 1 percent on Nr 6 j and less than 0.5 percent on total temperature. 



I.OxlO ' 2 



T, °R 


Figure 2k.- Thermal conductivity of chromel and alumel thermocouple wire and comparison 

with other "type K" wires. 


78 



APPENDIX B 


MEASUREMENTS OF HEMISPHERICAL EMISSIVITY FOR 
CHROMEL AND ALUMEL WIRES 

By William D. Harvey, Lemuel E. Forrest, 
and Frank L. Clark 

The use of thermocouple probes for surveying hypersonic flow fields over test 
models or tunnel wall boundary layers poses problems of compromise among such factors 
as ability to withstand high-temperature and high-velocity conditions, conduction and radi- 
ation errors, recovery characteristics, and often the time response of a probe. Usually, 
a thermocouple material can be selected and a probe designed and constructed to properly 
account for most of these factors. However, radiant -heat -transfer losses from a thermo- 
couple probe become especially important at elevated temperatures, and corrections due 
to radiation must be accurately calculated. 

Usually, sufficient data are available on the normal emissivity for individual base 
metal thermocouple wires at low temperatures (see, for example, refs. 60 and 61). How- 
ever, limited definitive data are available on total hemispherical emissivity for selected 
thermocouple wires at high temperatures with the exception of some results presented 
in references 61, 62, and 63, The indiscriminate substitution of the normal emissivity 
for total hemispherical emissivity in heat -transfer problems could lead to serious errors 
(paper no. 31, ref. 61). 

The purpose of this appendix is to present data on the total hemispherical emissivity 
for chromel and alumel thermocouple wires over a wide temperature range. These 
results have been utilized to calculate the radiation and conduction corrections by the 
method of appendix A for all data presented herein. 

The effective hemispherical emissivity e of an exposed thermocouple wire of 
length L and diameter d was experimentally determined by electrically heating the 
wire in an evacuated chamber and calculating the emissivity by use of equation (A5) 
or equation (A6) of appendix A. 


Apparatus and Procedure 

All chromel and alumel wires used in the present tests to determine the emissivity 
were 0.508 millimeter (0.020 in.) in diameter. In order to minimize conduction losses 
a large overall length -to -diameter ratio of L/d = 525 was used. A schematic of the 
wires and circuit used is shown in figure 25. A block diagram is included in this figure 
indicating the equipment used and method of data readout. Figure 25(a) shows the alumel 
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Figure 25.- Emissivity apparatus circuit. Dimensions are in inches (parenthetically in centimeters). 
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Figure 25*- Concluded. 


wire in place while figure 25(b) indicates the individual circuit (less the heating circuit 
and readout equipment) for the chromel and alumel wires. 

All junctions of the chromel -alumel wires were formed by resistance butt welding. 
There was no "bead" formed when the junction was butt welded. The test wire was freely 
suspended between supporting posts 0.27 meter (10.5 in.) apart. Chromel and alumel 
wires of 0.0762-millimeter (0.003-in.) diameter were attached 1.27 centimeters (0.5 in.) 
on either side of the test-section center wire to allow evaluation of heat conduction effects 
from equations (A5) to (A7) where = 0, dTj^/dt = 0, and = 0 for these tests. 

The results indicated that heat conduction effects could be neglected in the present emis- 
sivity measurements. 

The testing was done in a 0.61 -meter (2 -ft) square, steel, vacuum chamber which 
was evacuated to pressures of about 1.333 N/m^ (10 microns of mercury) maximum. 
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The test wire was electrically heated and measurements were taken of the voltage drop E 
across the central 2. 54 -centimeter (1-in.) test portion. A 60-cycle alternating current 
was used to heat the wire and this current I was measured with a calibrated ammeter. 

The thermal electromotive force of the junction at the center of the test section 
relative to the ambient temperature was measured with a manually balanced precision 
potentiometer. The potentiometer was isolated from the 60-cycle current by a low-pass 
filter. A back-to-back arrangement of electrolytic capacitors was used to prevent short 
circuiting of the thermal electromotive force by the secondary winding of the power trans- 
former. Effectiveness of the circuit was verified by observing that no direct-current 
potential difference existed across a shunt in series with the transformer windings. 

Tests and Results 

In tunnel tests, the survey temperature probes are exposed continuously to the test 
gas flow and over a period of time become "aged," that is, changes in the surface condi- 
tion of the thermocouple wires occur. It has been found and discussed in reference 58, 
for example, that both chromel and alumel wires are subject to oxidation in air above about 
1144 K (2060° R) but are still more resistant to oxidation than other commonly used con- 
ventional base -metal thermocouples. However, as found in the present test program to 
measure emissivity and as mentioned in reference 58, reducing conditions (that is, a test 
environment which has not been purged with high purity nitrogen prior to evacuation) have 
adverse effects and the chromel wire goes through a change in a surface condition called 
"green rot." Under this oxidizing condition, the chromel wire is oxidized to a much 
greater extent than the alumel wire. 

It was observed in the present tests that the chromel wire did darken with increasing 
exposure time and temperature level when the test chamber was not purged with a high 
purity nitrogen gas. The alumel wire did not darken as much as the chromel wire. The 
present results referred to as "oxidized" have been obtained by subjecting the wires to 
various indicated temperature levels for fairly long time intervals under a nonpurged 
vacuum system. A new wire could be kept "bright" by purging the vacuum chamber with 
high purity nitrogen gas prior to each test. The nitrogen gas used for purging contained 
less than 5 parts per million of oxygen and is the same gas used in the tunnel (ref. 30). 

The total hemispherical emissivity for bright and oxidized chromel and alumel wires 
is shown in figure 26 for successive runs on each wire. The present data are compared 
with available values for normal emissivity (ref. 64). The present results indicate 
increasing emissivity with temperature for both the chromel and alumel wires and the 
magnitude of emissivity is considerably greater for the oxidized wires than for the bright 
wires. The present data are generally somewhat different from previous values for nor- 
mal emissivity. 
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APPENDIX B - Concluded 


The results shown in figure 26 have been utilized in the data-reduction procedure 
to correct the total -temperature probe for radiation losses according to the method of 
appendix A. The difference (e^ - has been assumed invariant with wire exposure for 
the oxidized condition and the present results indicate that this assumption is valid. 

The dashed line curves in figure 26 were used for the values of required in 
the data-reduction procedure of appendix A. These curves are extrapolations to lower 
temperatures of the emissivity e at the edge of the boundary layer as determined 
from equation (A5) with the required condition that Tj = T^ e at y = 6. 
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